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ABSTRACT 


Under  this  effort  work  has  been  done  to  derive  and  use 
closed  form  expressions  for  generalized  impedance  parameters 
involving  rectangular,  conducting  patches  in  free  space.* 

A  Galerkin  model  for  the  interaction  between  two  wires,  a  wire 
and  a  patch  and  two  patches  has  been  developed.  Basis  functions 
sinusoidal  in  the  direction  of  the  current  and  triangular  in 
the  transverse  direction  have  been  used.  Farrar's  Integration 
technique  has  been  used  so  as  to  give  most  of  the  results  of 
the  model  in  closed  form. 


The  original  proposal  called  for  the  patches  to  be  parallel. 
This  restriction  has  not  been  applied  and  the  formulation  is 
general  with  respect  to  the  relative  orientations  of  the  trans¬ 
mitter  subsection  and  the  receiver  subsection. 
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CHAPTER  1:  FORMULATION  OF  THE  PROBLEM 


1.0  Introduction 


The  treatment  of  conducting  surfaces  using  the  method  of 
moments  has  been  a  complex  problem.  Wire  grid  models  have  been 
used  with  some  success  and  are  currently  available  in  the  GEMACS 
system  (Balestri  et.  al.  (1977)),  among  others.  Wire  gridding 
gives  good  results  for  far  field  problems  and  scattering  problems. 

If  impedance,  currents  or  near  fields  are  required,  or  if  the  surface 
contains  slots  comparable  to  the  size  of  the  grid,  the  results  from 
wire  gridding  are  poor.  Wire  grid  models  have  been  criticized 
because  of  questions  of  their  validity  (see,  for  example,  Lee  et. 
al,  1976) .  Surface  current  models  have  been  developed  by  Knepp 
and  Goldhirsh  (1972).,  Albertson  et.  al.  (1974),  Wang  (1974), 

Wilton  et.  al.  (1976),  and  Singh  (1977a).  Both  Knepp  and  Albert¬ 
son  et.  al.  use  the  magnetic  field  formulation  with  pulse  basis 
functions  and  point  matching.  Wilton,  et.  al.  use  the  electric 
field  integral  equation  with  pulse  basis  functions  and  point 
matching. 

Wang  uses  the  electric  field  formulation.  His  expansion 
functions  are  sinusoidal  in  the  direction  of  the  current  and  uniform 
in  the  transverse  direction.  The  same  applies  to  his  testing 
functions.  The  details  of  his  impedance  calculation  are  not 
given. 
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in  the  transverse  direction.  In  addition,  non-rectangular  shapes 
of  basis  function  patches  are  permitted. 

The  limiting  factor  in  all  techniques  mentioned  above  is  the 

large  amount  of  CPU  time  required  to  numerically  perform  up  to  four 

integrations  for  each  matrix  element  and  to  solve  the  resulting 

matrix  equations.  It  is  to  be  noted  that  if  £  wavelengths  is  some 

average  linear  dimension  of  a  surface  and  if  v  basis  functions  per 

wavelength  are  to  be  used,  the  core  required  for  the  matrix  is  of 
4 

the  order  of  (l\>)  .  The  CPU  time  required  to  fill  such  a  matrix 
is  also  of  the  order  of  (£v)^  and  the  CPU  time  required  to  solve 
the  equations  is  of  the  order  of  ( S. v )  ^ .  Such  high  order  dependence 
may  result  in  CPU  times  of  the  order  of  several  hours  for  problems 
of  moderate  size  (about  2  or  3  wavelengths  square) . 

Recently,  a  new  integration  technique  due  to  Farrar  (1978) 
et.  al.  has  been  published.  This  technique  has  been  used  to  perform 
some  of  the  integrations  in  closed  form  and  thereby  reduce  comput¬ 
ation  times  drastically. 

The  problem  has  been  formulated  in  the  next  three  subsections. 
Subsection  1.4  is  devoted  to  defining  the  notation  for  the  coord¬ 
inate  system  to  be  used  in  the  model  which  is  developed  in 
Chapter  2. 


•  i 


1-3 


1.1  Basis  Functions  for  Wires 

Consider  a  wire  subsection  i  extending  from  t  s  .1  to  l  =  The 

(current)  basis  function  $.(£)  is  given  in  terms  of  and  £2  as 


where 


= 


I" sin  k(£  ♦  £p 

sliTTcF^ 

sin  k(£2  -  *•) 
sin  k£„ 


4j'<  l  <  0 


0  <  £  <  £2 


elsewhere 


k  =  is  the  wave  number. 


(1.1.1) 


Lengths  £j,  £2  are  either  equal  or  one  of  them  is  zero.  Basis  functions 
4*^^ (£)  are  pictorially  depicted  in  figure  1.1.1. 

Figure  1.1.2  shows  a  wire  divided  into  several  basis  functions.  The 
presence  of  a  half  sinusoid  at  the  end,  if  that  end  forms  a  junction,  may  be 
noted.  The  half  sinusoid  implies  the  presence  of  a  point  charge  at  the  end. 
The  expressions  in  the  next  sections  take  this  point  charge  into  account. 

<(>(£)  is  dimensionless.  The  current  on  the  wire  segments  in  the  configuration 
is  given  by 

i  l  i  x,i 

where  1^  is  the  current  on  wire  segment  i,  u^  is  the  unit  vector  in  the  £ 

i 

direction  for  wire  segment  i. 


1  / 
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1.2  Basis  Functions  for  Surfaces 

Consider  a  surface  subsection  i  carrying  current  in  the  u  direction  and 
extending  y  =  -£3  to  £^  and  z  =  -£^  to  The  current  density  basis  function 

for  J^ly.z)  is  given  in  terms  of  £^,  £3  and  as  iJk Cy) (Z)\LZ  where 


£3  ♦  y 
£ 


£3  <_  y  _<  0 


I'ify) 


0  <  y  <  £„ 

—  4 


(1.2.1) 


0  elsewhere 

Lengths  £3  and  £^  are  either  equal  or  one  of  them  is  zero.  It  is  to  be 
noted  that  the  (x,y, z)  coordinate  system  is  local  to  the  surface  subsection 
under  consideration  and  is  defined  in  such  a  way  as  to  have  the  current  in  the 
subsection  flowing  in  the  x  direction,  to  overlapping  subsection  carrying 
current  in  the  orthogonal  direction  is  a  part  of  the  set  of  basis  functions. 

Thus  the  method  of  moments  solution  may  be  expected  to  contain  current  density 
components  in  any  direction. 

Depending  upon  the  values  of  £^,  ^3  anti  nine  different  types  of 

basis  functions  are  possible.  These  and  their  two  dimensional  representations 
are  shown  in  figure  1.2.1.  Figure  1.2.2  shows  the  complete  set  of  basis 
functions  for  a  surface  ABCD  whose  edges  AC  and  CD  are  connected  to  other  planes. 
The  current  density  on  the  surface  segments  in  the  configuration  is  given  by 


*ln  this  report,  the  term  "surface  patch"  implies  a  surface  expansion  function 
and  not  a  physical  section  of  the  surface.  It  is  assumed  that  a  local  coord- 
ate  system  has  been  arranged  so  as  to. have  the  patch  in  the  yz  plane,  the 
z  direction  being  the  direction  of  the  current. 


/ 


—Vjtrf 1 


Figure  1.2.1:  Pictorial  (3D)  representation  of  all  possible  current 
density  basis  function  types  possible.  Some  of  the  corresponding  2D 
representations  are  drawn  below  the  3D  representation  and  some  others 
are  drawn  above.  The  2D  representations  and  the  corresponding  3D 
representations  are  connected  by  dotted  lines. 


the  orthogonal  current  densities  separately. 
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£  J  =  E  <My)<Mz)J.u  (1.2.2) 

i  x  i  1  A  A  zi 

where  is  the  contribution  of  the  current  density  basis  function 

ip. (y)4>^(2) u  to  the  total  current  distribution.  All  J.  are  unknowns  in  the 

*  _  i 

problem  and  have  units  of  amps  m  .  and  are  dimensionless  functions. 


•  / 


—  i 


1-10 


1.3  Method  of  Moments  Formulation 

The  geometry  of  the  problem  under  consideration  may  involve  wires  and 
surfaces.  The  unknown  in  the  case  of  wires  is  current  whereas  the  unknown 
in  the  case  of  surfaces  is  current  density.  The  method  of  moments  yields  the 
following  system  of  linear  equations: 


[ZWW]  [ZWS] 
[ZSW]  [ZSS] 


[I] 

[VW] 

_[J]  - 

.  1  Vs] 

(1.3.1) 


The  superscripts  W  or  S  indicate  whether  the  submatrices  in  question 
pertain  to  wires  or  surfaces.  Furthermore, 

vi  =  / EflO.u^  <Di(2-Dd2.  (1.3.2) 

and  * 

visff  £(y.z)-uz  4>i(z)ij)i{y)dzdy  (1.3.3) 

C 

where  E(.)  is  the  incident  electric  field.  It  is  to  be  noted  that 

ww 

has  units  of  volt-m.  It  may  be  noted  also  that  whereas  Z  has  units  of  ohms, 

WS  SW  SS  2 

Z  and  Z  have  units  of  ohm-m  and  Z  has  units  of  ohm-m  .  The  mutual  impe- 

WW 

dance  Z^.  between  two  wire  segments  (at  y'=0  and  at  n=0)  shown  in  figure  1.3.1 
is  given  by 


yfyj  r  2  r  2  _ 

Zij  =  J  c J  fG(0,c;0,z').  ^(z')*^,  dz'dc 


a.  3. 4) 


—  s. 
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Here,  G(0,j;;0,z')  is  the  free  space  dyadic  Green's  function  and 
-2 

units  of  ohm-m 
a  surface  is  given  by 


has  units  of  ohm-m  2 .  The  mutual  impedance  Z^  between  a  wire  and 


-a 

j 


zijaJ  /  ,  G(0,c;y*,z').^j(z')uzf  dz'dy'dr, 

(1.3.5) 


Similarly, 


2SW  4  ^  (n)  /"  2  <J>.(t)u  .  /"  2  G(n,C,’0,z')  .(}> .  (z')u  ,  dz'd^dn 

(1.3.6) 


and  finally, 
rZ 


hu  Ve.  ■  •'-««  3  -V 


G  (n,:,;y’  ,z' ' 


#  j  (z ’ ) uz , dz ‘dy 'd^dh 


'I.3.' 


In  the  above  equations,  the  element  i  has  been  used  as  the 
receptor  element  and  j  has  been  assumed  to  be  the  source  element. 
Because  of  the  symmetry  of  the  equations,  this  choice  is  arbitrary. 
In  the  next  chapter,  techniques  for  evaluating  the  impedance  quant¬ 
ities  are  studied. 


>  / 
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1.4  A  Notation  for  the  Coordinate  Systems 

Figure  1.3.1  shows  a  receiver  and  a  transmitter  for  which 
the  impedance  model  is  to  be  developed. 

In  general,  the  coordinates  and  distances  from  a  point 
(0,y',z')  on  the(x,y,z)  system  to  a  point  (0,n»C)  on  the 
(x,n*c)  system  are  of  interest.  Of  interest  also  are  some  of 
these  values  when  one  or  toore  elements  of  the  set  <y*n»C>  are 
zero.  A  three  subscript  notation  is  used,  each  subscript  being  a  0 
or  a  1.  If  a  subscript  is  0,  the  corresponding  variable  in  the  set 
is  0.  If  the  suoscript  is  1,  the  corresponding  variable 
in  the  set  is  present.  Thus,  rQ00  is  the  distance  between  (0,0,z') 
on  the  (x,y,z)  system  and  (0,0,0)  on  the  (x,n,C)  system.  Simil¬ 
arly,  is  the  distance  between  (0,0,z')  on  the  (x,y,z)  system 

and  (0,n,S)  on  the  (x,n»C)  system.  This  notation  is  extended  to 
allow  the  value  (-1)  for  the  subscripts.  Under  this  extension, 
r01(-l)  the  ^stance  between  (0,0,z‘)  on  the(x,y,z)  system 
and  (0,  n,  -j;)  on  the  system.  Thus,  when  the  subscript  is 

I 

-1,  the  corresponding  variable  in  the  set  <y  ,n,C>  is  negated. 

The  distance  from  (0,y',z')  to  (0,n,C)  is  given  by  the  vector 


-  / 


.  i. 
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-111 


where 


=  (X111'  *111'  *111>  =  £000  +  V  y  +  -nn  x  V 

(1.4.1) 


U  1 

-y 

=  (0, 

1 

H-* 

O 

u 

=  (a  , 

a 

.  a 

-n 

xn 

yn- 

and 

Hr  = 

(a 

'  V 

(1.4.2) 


In  the  above  equation,  Eqqo  ^-s  the  vector  from  the  point 
( 0 , 0 , z 1 )  on  the  (x,y,z)  system  to  the  origin  of  the  (x,n,U  coord¬ 
inate  system.  a  ,  a  ,  a  ,  a  ,a  ,  a  are  the  directional 

y  s»  zc,  xrj  yr)  zn 

cosines  between  the  respective  coordinate  axes.  Following  the 
properties  of  directional  cosines, 


+  a 


yc 


+  a 


z? 


1, 


(1.4.3a) 


+  a 


yn 


+  a 


zn 


1. 


(1.4.3b) 


Quite  frequently,  it  is  desired  that  some  of  the  distances  be  expressed 
as  a  quadratic  trinomial  in  one  of  the  variables.  A  notation  for  syst¬ 
ematically  doing  so  is  presented  here.  The  starting  point  under  this 
notation  is  always  r^^  which  may  be  written  as 


1  / 
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rm  'n'c>  =  R  =  (^no  +  V>-  (r110  +  U rp 
=  (r^1Q  +  ii0 ' - +  l2)  =  (Rg  +  R1  C  +  ^2> •  (1.4.4a) 

The  arguments  of  r...  may  be  omitted  if  they  are  the  variables 
y',n  and  5.  The  £  in  equation  1.4.4  is  a  superscript,  not  power, 
Similarly, 


rlll  ~  (R 


•0  +  R1  y'+  y'2>  and 


(1.4.4b) 


rin  =  (Ro  +  Rin  +  n2) 


(1.4.4c) 


Each  of  the  coefficients  of  the  quadratic  trinomials  Ri ( ’ *  can 
themselves  be  written  as  polynomials  of  degree  (2-i)  in  other  var¬ 
iables.  Thus,  for  example, 

rll0  (y’'R)  k  rq  =  (£100  +  VK  (Sioo  +  V°  = 


(rlOQ  +  ^IQO'V  n  +  n“>  4  Rnn  +  H  +  H 


=  “00  01 


(1.4.5) 


Table  1.4.1  shows  all  such  expansions  that  have  been  used  in 
this  report.  They  have  been  derived  using  a  development  similar  to 
1.4.4a,  above. 

Whenever  the  ensuing  text  calls  for  a  coefficient  one  of 
whose  subscripts  is  2,  the  coefficient  is  assumed  to  have  the  value 
1.  Whenever  the  sum  of  the  subscripts  is  greater  than  2,  the  value 
of  the  coefficient  is  0. 

The  vector  £^11  =  ^xlll'  ^111'  ^  also  importance  in  the 
calculations  of  impedance  in  cases  where  a„_  5*  1. 


■-fW+'T' 


t>»»» 
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Pill  =  £000  +  V  y'  +V  +  V 


(1.4.6) 


where 


V  =  (0'  -1'0)'  -n  =  (axn'  ayn '  0) '  ^  =  (axe'  V'  0) * 

(1.4.7) 


Expansions  for  p  that  are  similar  to  those  given  above  for  r 
are  possible  and  are  represented  by 


111 


»m  ’  A  *  <Ao‘  +  *i<  *  <2>  i 


(1.4.8) 


The  expansions  are  similar  to  those  for  r^^  except  that  none 

of  the  coefficients  of  the  trinomial  is  known  apriori  to  be  1.  It 

2 

is  convenient  to  express  the  trinomial  using  the  multiplier  v  (as 
has  been  done  in  equation  1.4.8),  when  expansions  about  C  are  involved. 
It  is  more  convenient  to  express  it  using  three  coefficients  (as  has 
been  done  in  the  expansions  for  Ag  in  table  1.4.2)  when  expansions 
about  n  or  y‘  are  involved.  All  expansions  are  listed  in 
table  1.4*2.  It  is  noted  that  most  expansions  do  not  exist  when 


v  =0 . 


Finally,  the  quantities  a  and  cos  0^.^  are  exPressed  as  a 


function  of  e  thus: 


PC 


_  £nr  -e 
-£111! 


(£ll0  +  v,C) .uc 


Be 


(aJ  +  aJc  +c2)1/2‘ 


vc  (i  Ai +  c) 


2^*5 


Y5  I  (Aq  +  Aip;+  c2) 


?cl 

(A 


n  +  AiC+  C2)** 


(1.4.9) 


cos  (0.^^)  =  HI 


’111 


g110  +  azcC 

rlll 


(1.4.10) 
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TABLE  1.4.1:  Expansion  of  the  Distance  Parameters  into  Quadratic 
Trinomials 


Coefficient  Expansion  Interpretation 


R 

Rg  +  R^  ^  +r>2 

2 

rm 

Ro 

Rcy’  +  Rcy'  v.  +  v, 
R00  +  R01  y  +  y 

2 

2 

r110 

RC 

Ro 

Roo  Roi  n  n 

2 

r110 

R1 

Rsy’  +  R?y*  v. 

R10  +  R11  y 

2-110  -r. 

R^^  +  R^^  n 

R10  Rn  n 

2— 1 10 "  —  r. 

K00 

RCy,n  cy'n  + 

Kooo  001  n 

2 

n 

2 

r010 

R^y’ 

R01 

Rcy'n  .  Rcy'nn 
Roio  +  Ron  n 

2  —0 10*  — y 

R^n 

R00 

Rcny*  +  Rcny'  v.  + 
Rooo  +  Rooi  y 

.2 

y 

2 

r100 

RS0 

R01 

Rcny'  .  Rcny'  , 
R010  +  R01 1  y 

2-100 " -n 

R^y’ 

K10 

R^y'n  +  R^y'n  n 
R100  R101  n 

2— 010 "  —  r. 

R^y’ 

R11 

Rcy'n 

Ril0 

2ay '  r. 

I 

•  / 
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TABLE  1.4.1  (continued): 


Coefficient 


Expansion 


Interpretation 


R 


R 


R 


R 


R 


R 


R 


R 


R 


R 


R 


cn  Rcny'  .  R5ny’ 

10  100  R101  y 

2-100*-f; 

sn  RCny’ 

11  110 

2V  <-«l 

?y’n 

000 

rooo 

cy'n 

001 

2-ooo  *-n 

cy'n 

010 

2— 000"  -y  ’ 

cy'  n 

on 

2a  , 
ny ' 

Cny ' 

000 

2 

rooo 

cny ' 

001 

2— 000 y 1 

cny' 

010 

2-ooo ■ — n 

?ny ' 

Oil 

2a  , 
ny 

5y'n 

100 

2-000‘-? 

11 


2  *■■*■•)*■* 


_  t 
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TABLE  1.4.1  (continued): 


Coefficient 


Expansion 


K 


ty'n 

101 


R 


cy'n 

no 


,cny* 

100 


R 


cny‘ 

101 


Interpretation 


2V0) 


2a  ,  „ 

y'c 


2^ooo 


2a  , 

y'c 


2an5,-0) 


I 
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TABLE  1.4.2:  Expansion  of  p  in  Quadratic  Trinomials 


Coefficient  Expansion 

Interpretation 

A 

vc  (Ao  +  Aic  +  ;2) 

2 

plll 

4 

Acy*  .  Asy’  v.  .  Asy’  v.2 

A00  01  y  02  y 

P110/v2 

*5 

Ai’r'  +  A^n  +  A^n* 

00  01n  02n 

P110/v2 

Aey'  +  A?y'  . 

Aio  An  y 

2^  no-Vvc 

A?n  +  A?nn 

A10  Alln 

2£  uo-V’! 

41' 

Acy'n  +  Asy'n  +  Acy'n  * 

000  001  n  002  n 

p010/vc 

A5y' 

A01 

Asy'n  +  Acy'n 

A010  A011  n 

2£„.  .h-'t 

A^y' 

A02 

Acy’n 

020 

v2,/v2 

y’  c 

4S 

Acny'  .  Acny '  v,  .  Acny'  ,2 

A000  A001  y  A002  y 

2  2 
piooA 

*sr 

AW  +  Acny’  , 

A010  on  y 

2£  100 ‘-n  /vr 

A02 

Aw 

020 

v2/v2 

n  c 
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Table  1.4.2  (continued): 

► 


Coefficient 


Expansion 


Interpretation 


ACy' 

A10 

yv^y'n  +  A?y’nn 

100  101  n 

^oio^c 

ACy' 

11 

AYy'o 

110 

2V*V\ 

A^n 

A10 

ASny'  .cny' 

100  101  y 

2Lm’Uv\ 

A11 

ACny' 

A110 

2Y.n-Y.^/\(=  °) 

Acy’n 

000 

pooo/vc 

A^v'n 

A001 

2£ooo'Vvc 

A^y'n 

002 

v2/v2 
ir  e 

A?y,ri 

010 

2£000’-y,/v;2 

A^y'n 

on 

2vn.^,/v2 

A^y'n 

020 

V/v? 

A^y 

000 

P000/v' 

:*S"**t* 
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TABLE  1,4.2  (Continued): 


Coefficient 

Expansion 

Interpretation 

Asny' 

001 

^OOV7^  ' 

Acny' 

002 

2.  2 
vy ' /vc 

Acny* 

A010 

2tO00  %/v\ 

A^ny' 

Oil 

2vy, .vn/v^ 

ACr,y' 

A020 

2/  2 

Vvc 

Asy,n 

A100 

2M)00*V  \ 

A^y,n 

A101 

2VVvc 

AS’n 

110 

2— y  '  *y5/v^ 

Acny’ 

100 

2t000*-/v/ 

A^y’ 

101 

2yy .  -y r/v* 

A^ny' 

110 

2-r\'-r/\ 
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CHAPTER  2.  CLOSED  FORM  SOLUTIONS  FOR  IMPEDANCES 


2. 1  Coupling  Between  Two  Surface  Segments 


The  mutual  impedance  between  two  surface  segments  i  and  j  is 


given  by  (see  equation  1.3.7) 
i. 


zf®  =  f  0H(n)dn  f  0.  U)u  .  d  tf  0s(y’)dy'  f  G(n,C;y',z' 


<J)  (z‘  )u  ,dz' 
J  Z 


(2.1.1) 


To  calculate  this  mutual  impedance  we  first  consider  the  mutual 
impedance  between  the  surface  segment  i  and  a  wire  on  segment  j  at  a 
given  value  of  y' .  This  mutual  impedance  is  given  by 


/*■  4  *2  *2 

0,  (n)dn  f  0  •  ( c) u  »d£ f  G(n,t;y',z' 
■8.-1  1  J-SL,  1  J-SL! 


.2  £2 

i 

'1  ~1 
From  2.1.1  and  2.1.2  it  is  easy  to  see  that 

Z?j  =  f*  Vy,)zij<y',dy'- 

Similarly, 

ZSS  =  f  4  0i(n)z”j  (n)dn 


)  .0  •  (Z  ' )  U  ,dz ' 
J  z 

(2.1.2) 


(2.1.3) 


(2.1.4) 


ws 

where  (n)  is  the  mutual  impedance  between  a  wire  on  surface 
segment  i  at  a  given  value  of  n  and  the  surface  segment  j.  Equa- 

ss  ws 

tions  2.1.3  and  2.1,4  portray  the  hierarchy  between  Z  and  Z 
sw 

and  Z  in  that  the  integration  of  the  latter  with  respect  to  one 


•  / 
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ws  sw 

of  the  coordinates  gives  the  former.  Thus,  finding  Z  or  Z  is 

ss  WS  sw 

a  step  towards  finding  Z  .  The  computation  of  Z  and  Z  is 

SS 

treated  in  the  next  section.  The  problem  of  evaluating  Z  using 

equation  2.1.4  is  revisited  in  Section  2.5  following  the  computati 
wc:  cw 

of  Z  and  Z°  in  Section  2.4. 


>  l 
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2 . 2  Coupling  Between  Wire  and  Surface  Segments 


Using  arguments  similar  to  those  used  in  Section  2.1,  it  is 


easily  seen  that 


*'4 

:sw  •  /, 

ID  ID 


(2.2.1) 


ww 

where  Z^(n)  is  the  mutual  impedance  between  a  wire  on  segment  i 

sw 

at  n  and  the  wire  j.  Z^j(y’)  defined  in  equation  2.1.2' is  given 
by 


z^(y')  =  f  iK(n)  z^(n,y’)dn 

A J  •/  _  p  A  x  J 


(2.2.2) 


IVW 

where  ij  Cn,y')  ttie  mutua^  impedance  between  a  wire  on  segment  l 
at  n  and  a  wire  on  segment  j  at  y' . 


Similarly, 

zii  =  f 4  'My,)zii(y,)dY' 

1 J  J  _  p  »  J  1 J 


( 2 . 2  *  3a  ) 


•3"»  ■ 


0  • 

*4 


r 


.WW  , 


(y  * ) z .  .  (y’ , n) dy 1 . 

o’  J  !3 

*3 


(2.2.3b  ) 


'  /  _ -  *•«**»»■»  «ai 
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WS 

Equations  2.2.3  and  2.2.2  portray  the  hierarchy  between  Z 
sw  ww 

and  Z  and  Z  ,  in  the  sense  that  integration  of  the  latter  with 

respect  to  y'  or  n  gives  Z^  or  Z^W  respectively.  Thus,  finding 

ww  WS  SW  WW 

Z  is  a  step  towards  finding  Z  or  Z  The  computation  of  Z 

WS 

is  treated  in  the  next  section.  The  problem  of  evaluating  Z  or 
SW  WW 

Z  from  Z  is  revisited  in  section  2.4. 


*  *  -'mfi 


Equations  1.3.4  through  1.3.7  contain  the  term 


E, 


(i)/Jj 


(.).z').dj (zl)uz.dz'  2.3.1 


as  their  innermost  integral.  This  quantity  is  the  r  directed 
component  of  the  electric  field  at  (0,n,i‘.)  due  to  a  unit  current 
on  subsection  j  and  has  units  of  volts  m  3  amp-^.  The  result  of 
the  integration  is  a  vector  whose  dot  product  with  u  gives 
E .  ( i ) / 1 j .  A  closed  form  solution  has  been  derived  by  Schelkunoff 
and  Friis  (1952)  and  restated  in  convenient  terms  by  Singh*  (1977) 
When  I.;  is  1,  E.  ^  is  given  by 

E-.  (i)  (*  -  *  )  =  -  TF.  +  az,-*z  +  '-2*3-2' 

where 


V  (v '  ,t|,r.)  =  - 


4  ii  io  i 


l’  "  =-  j  ,  0  ,  V  2 


,-jkr 

,  Is'  (z*)l  -r— 


111 


(2.3.3) 


(.  .  .  ) 


_ L 


4  ii  lot 


E  , 


«  '=“^,0,  V2 


h<z'  )i  o:3!fiii  r  ,  ^ 

•  L  J  ^ui  .  ' 


♦Pages  40-42. 
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and 


E  ( . . . ) 

p 


4nujep 


(...) 


1 


,0,* 


A 

£* 


(z1)  ]exp(-jkr11;L) 


+ 


.  3 

4ttu)£  p 


(...) 


A 


cose 


in 


exp(-jkrni) 


(2.3.5) 


where 

A  (.)  =  (.)  -  (.)  (2.3.6) 

t '  i'+O  i,*-0 


Since 

I j  = 1  and  is  given  by 


ww 

Z  is  independent  of  I . , 
13  3 


it  is  computed  using 


.WW 

i  j 


(y* ,n) 


M 

j  (y’,n,e)  dc 


(2.3. 7) 


WW 

Since  the  arguments  of  the  functions  Z  ( . . )  ,  <(.(.)  and 
E  ( . . . )  provide  the  context  for  the  subscripts  i  and  j,  these 
subscripts  are  dropoed  in  the  ensuing  discussion. 

Substituting  2.3.3,  2.3.4,  and  2.3.5  into  2.3.2  and  then 
substituting  2.3.2  into  2.3.7,  the  following  equation  is  obtained: 


ZWW 


(y' /n) 


4  WW  WW 
Z  b.  Zt  (£' ,y'  ,n) 
t=l  c 


>  i 


(2.3.8) 
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where 


iir  J, 


3<fr(z’)1  . 

^  _  i  I  » 


b""  =  _ 2_  A  2±±±_t 

2  4ttu)£  P,  3  z 1 


.WW  k  . 
b3  =  w  J, 


uWW  _  j  A  3»  (Z1  ) 
4  ~  4 it (ue  P,  3z  * 


(2.3.9a) 


(2.3.9b) 


(2.3.9c) 


(2.3. 9d) 


and  where 


zV”  U\y\n)  = 


_WW  , - ,  ,  , 

z2  U  »y  »n)  = 


zT  (t'.y'.n)  = 


„ww  .  „  ,  .  . 

z4  U  ,y’  ,n)  = 


3  /exp(-j 


~  jkr,  , ,  )> 


exp(-jkr, . . ) 

a  He)  - r - dc 

ZC  rlll 


'2  exp (-jkr. . . ) 

a  He)  - dc 

PC  plll 


(2.3.10a) 


(2.3.10b) 


(2.3.10c) 


cose  exp(-jkrin) 

a  He) - — - —  d5  (2.3.  lOd) 

pe  p1u 
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In  the  above  equations,  Z^  is  the  contribution  of  the 

ww 

current  discontinuity  on  the  source  wire.  Z^  is  the  contribution 

of  the  z  directed  field  generated  by  the  source  wire.  This 

contribution  is  zero  when  the  two  wires  are  orthogonal  to  each 
WW 

other.  Z ^  is  the  contribution  of  the  p  directed  field 

generated  by  the  source  wire  segment  when  it  has  a  current 

discontinuity  at  one  end.  The  contribution  of  Z^  is  zero 

whenever  the  two  wire  segments  are  parallel  or  when  the  source 

WW 

segment  has  no  current  discontinuities  on  it.  Z^  is  the 
contribution  of  the  p  directed  field  generated  by  the 
source  wire.  This  contribution  is  zero  when  the  two  wires  are 
parallel. 


Theorem  2.3.1 

Each  of  the  component  impedances  Z^  through  Z^  can  be 
written  as  linear  combinations  of  the  impedance  quantity 
( 5. 1  ,y'  ,(i,m,s,b)  where 

exp(-iVr  ^  dc 

(2.3.11) 


ri-. 


)  = 


<P  (?)  (Ac 


C 

+  AiC 


+  C2) 


*4V 


111 


1  / 
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Proof 


The  proof  constitutes  finding  the  linear  combinations  of 
. )  that  give  the  component  impedance  in  question. 


zf(...)  = 


M 


aT 


-l: 


exp(-jkrin) 


111 


dC 


<p  (?)  (l  +  jkrni)  (is  +  0  r”^  exp(-jkri;L1) 


11  1-h  WW 

Z  Z  (jk)9  Or’  )  zn  ( .  .  .  0  ,h,g-3)  (2. 

g=0  h=0  0 


_ww .  .  _ww  ,  .  .  , . 

Zj  (...)-  azc  z0  (.. .0,0,-l) 


(2. 


-*i 


'2  0  Aj  +  ?) 

(Ao  +  a?  S  +  C2) 


4>  ( C )  exp  (- jkrxll)  d? 


2  [h  *\  )  z  (..  .  (-2) ,h,  0) 
h=0  0 


(2 


-  /  ' 


d? 


3.12a) 


3.12b) 


.  3.12c) 
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rww 

'4 


.) 


0  A i  +  c) 

(Ao  +A^C  +  s2) 


<t>  (c) 


^z110 

rlll 


exp(-jkr1:L1)  dc 


1 

I 

h=0 


1-h 

) 


Z'l  (  ...  (-2)  ,  (q  +  h)  ,-l) 


(2.3. 12d) 


End  of  Proof 

The  use  of  theorem  2.3.1  suggests  figure  2.3.1  as  a  flowchart 

WW 

of  a  procedure  to  calculate  Z  .  Note  that  each  of  the  4  blocks  on 

WW 

the  right  hand  side  calls  for  ZQ  . 

WW 

A  closed  form  solution  for  ZQ  is  developed  next.  In  developing 
this  solution,  it  is  assumed  that 


,3*  - 


5 

i£o 


(jx)1 


(2.3.13) 


where  x  is  less  than  tt/2.  According  to  Abramowitz  and  Stegui  *  . 
this  approximation  will  result  in  less  than  .09  percent  error. 
More  terms  may  be  included  for  greater  accuracy.  However,  the 
errors  developed  during  solving  the  method  of  moments  equations 
are  likely  to  undo  any  potential  benefits  of  greater  accuracy. 


*Page  76: 

aQ  =  1  a1  =  1  a2  =  .  49670  a3  =  .  16605  a4  =  .03705  a&  =  .00761 


'v. 
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The  above  integration  range  is  divided  into  the  (-ij  to  0) 
and  (0  to  i2)  ranges.  In  the  first,  a  change  of  variables  from 
C  to  -z.  gives 


z  r  a  sin(kU  -c) )  .  .  r 


X  C 


(rll(-l fj  exP(^k  rll(-l)0)  dc 


(2.3.14) 


2  (-1) CTS  exp(-jkr. lft) 


=  r. 


no' 


=1  2 j  <?in(ki0) 


(*0  +  <-l>°  A$  C  +  C2)  '  CS  r“1(_1} 


exp 


jk  (ia  -  c  -  r 


ll(-l)a  +r110 


jj-expfjk^  +  s-r^.^+r^ 


dc 


(2.3.15) 


The  argument  of  the  exponent  is  less  than  i a  in  the  domain 
of  integration.  Application  of  equation  2.3.13  gives 


•  / 


Calculate 

distances,  currents, 
current  slope, 
relative  orienta¬ 
tion  of  the  wires 


/  Current  \ 
discontinuity 
on  the  wire? 


Current 

slope  continuous 
or 

.^perpendicular 

wires 


Current  ^ 
discontinuity 

or 

non-parallel 
wires  / 


Current  \ 
slope  continuous 
or 

"-vparallel  wires/ 


equation  2.3.8 


Figure  2.3.1:  A  Procedure  for  Calculating 
7WW 
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2  (-1) JS  exp (- jkr  . n)  5  ■ 

•■>  ■  1  — 5T -.t.,— « -  ,Lai«k' 


m, 

zo  1  •  •  ‘ 


0  =  1 


i=0 


(Afl  +(-l)a  A  i  ;  +  C2)i5V  CS  ,  ,  x  a 


11  (-1) 


( S.<j  ~  C  “  r 


ll(-l)  +  r110 


)i  -  (~«a  +  A  "  rll(-l)°  +  rllo)  1 


dC 

(2.3.16) 


which  by  repeated  use  of  the  binomial  theorem  gives 


2  (-l)as  exp(-jkrno)  5  .  i  m  i-m  ( i 

)-  E  - yy-TV. . C-. -  s  a^jk)1  E  E  E  [m 

i=0  m=0  n=0  p=0 


Zm( 

zo  ( 


2j  sin  kla 


i-m 

P 


x  (-l)m  fl-(-l)i4n+P 

. 


.  l-m-p  p 

l  r  r  “ 


110 


f  0 

j  ( A o  +  (-l)a  A i  C+C2): 


v  s  _a  a  rm-n  r“  dr 

c  rll(-l)  C  rll(-l)CT 


(2.3.17) 


The  quantity  inside  the  integral  is  normalized  (i.e.,  the  inte¬ 
gration  limits  are  made  to  be  0  and  1)  by  writing 
C  =  la  x  and  noting  that  =  (Rg  +  C  +  C2)5*  to  give 


— ■  t 
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zww 

zo  (-** 


•  • )  ~  r 

0=1 


2  < -l)os  exp(-jkr,10)  *®+1  (a0?) 


%v 


2j  sin(JtA0)r 


110 


r  a,  Ok*,)1  £  ;  ,  fj]  M  (*-)  <-i>»  fi-^,i+n+P] 

i=0  1  m=0  n=0  p-0  W  t”M  p  J  (  J 


frllo')n+p 
|~J  ZA2v(s+m-n)  (n+a)  ^Aj ' 


(2.3.18) 


where  the  elements  A^  of  the  vector  £a}  and  R^  of  the  vector  Ol 
are  given,  for  0«is2,  by 


Ai  -  h0  ('-110  i0)l#i 


’  *1  * 


(2.3.19) 


and  where 

ZA2NPQ 

where 


A 

<[A],[R])  =  f  (A  (x)  )  **N  x  P  (R(x))*°dx 
J  n 


(2.3.20) 


A(x)^  An  +  Ai  x  +  Aox  /•  R(x)§  RQ  +  R1  x  +  R0x  .  When  N=0, 


-  0  1 
2...  reduces  to  the  simpler 

2  f1  1 

1  r»A  / r  nl  l  _  I  * 

~  J  n  x 


2A1PQ  ([r]) 


(R  (x)  )  ^  dx  . 


(2.3.21) 


_  1 
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Techniques  for  evaluating  in  closed  form  are  detailed  in 

Appendix  Al,  techniques  for  evaluating  Za2NPq  are  detailed  in 
Appendix  A2. 

A  few  comments  regarding  terminology  are  in  order;  The  unsuper- 
scripted  and  ,  being  elements  of  the  vectors  [A]  and  [R]  are 
distinct  from  the  superscripted  variables  defined  in  Section  1.4 
(some  of  the  unsuperscripted  A^  and  R^  are  related  to  the  superscripted 
variables  of  Section  1.4  by  equations  2.3.19).  Furthermore,  the 
symbols  [A]  and  Qr],  being  aggregates  of  the  quantities  related 
with,  respectively,  the  two  dimensional  distance  (p )  and  the  three 
dimensional  distance  (r)  between  the  receiver  and  the  transmitter/ 
are  vectors  when  the  wire-to-wire  interaction  is  under  investigation. 
These  symbols  represent  two  dimensional  matrices  when  wire-to- 
surface  and  surf ace-to-wire  interactions  are  considered  and  three 
dimensional  matrices  when  surface-to-surface  interactions  are 
considered. 

In  summary,  the  mutual  impedance  between  two  wire  segments 
with  sinusoidal  basis  functions  is  composed  of  four  impedance 
quantities  (equation  2.3.8)  each  of  which  can  be  written  as  linear 

ww 

combinations  of  the  impedance  quantity  (theorem  2.3.1, 

WW 

equations  2.3.12).  The  impedance  quantity  can,  in  turn, 

be  written  as  a  linear  combination  of  the  quantity  2A2NPQ  for 
various  values  of  N,  P  and  Q  (equation  2.3.18).  Whereas  techniques 
for  evaluating  AA2NPQ  are  exan,i-ned  in  Appendix  A2,  it  is  noted  here 


/ 
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that  za2Npq  is  a  dimensionless  quantity.  Therefore,  (from  equat¬ 
ion  2.3.18),  Z^has  dimensions  of  (length)  s+v+a+1  which  is 
dimensionally  consistent  with  equations  2.3.8  and  2.3.12.  Further 
details  regarding  the  computation  of  wire-to-wire  impedance  are 
treated  in  Appendix  D. 
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2 . 4  Coupling  Between  a  Wire  Segment  and  a  Surface  Segment 
Revisited 


From  section  2-2, 


ZWS ( n )  « 


*< y')ZWW(y',n  )  dy'. 


(2. 


zSW(y')  - 


ww 

'('(t)z  (n/y1)  dri 


(2. 


-I- 


Replacing  the  dummy  variable  y1  or  n  by  -y‘  or 
the  negative  part  of  the  integration  domain. 


-n  in 


l ' 


zWS ( n )  =  t 

t  =  3 


T  2™  ( (-1)  Ty 1  ,n)  dy’ 


i  ’ 

x 


(2. 


and 


^SW.  .X 
Z  (y  ’ ) 


4 

Z 

C  =  3 


Zm  (y’  ,  (-l)£n)  dn 


(2. 


Using  the  substitutions  y*  = y  and  n  =  n^y  respectively, 

the  above  equations  may  be  rewritten  as 

1 


ZWS(n)  =  z  Z  (-l)e  l'  [  ye  Z^((-l)Tt;y,n)  dy 
t=3  e=0  T  ' 


(2. 


and 


1  / 


4.1a) 


4.1b) 


4.2a) 


4.2b) 


4.3a) 
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± 

ZSW(y')=  Z  E  (-l)f*f  [  yf  Z^y'.t-l)5  ft  y)  dy 
r  =  3  f=0  5  J  4 


(2.4.3b) 


Here  the  variable  y  is  a  dummy  variable  of  integration  and 
is  unrelated  to  the  y  used  in  the  coordinate  systems. 

Equation  2.3.8  can  be  substituted  into  equations  2.4.3  to 


give 

ws  4  1 

ZWS ( n )  =  E  E  (-1 )~i 


t  =  3  e=0  T  ,0,J.’2  t=l 


E  b™  Z^  n) 


(2.4.4a) 


and 


ZSW(y,)  •  i  i  (-1)1  S  (-'.yM-ll'tJ 

£  =  3  f =0  ?  l  ,0 ,  S,’2  t=l  tf  S 


4  1 


4  _WW  „sw 


(2.4.4b) 


whe  re 


z ™  n)  =  j  ye  Z™  (ft',ft;y.n)  dy 


and 


z^  (t'»y'#t?)  =  j  yf  zf  (i’,yMcy)  dy 
o 


(2.4.5a) 


(2.4.5b) 
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It  is  to  be  noted  that  equations  2.4.5  represent  definitions 

WS  SW  r 

of  the  quantities  Zfcg  and  Zfcf  and  that  the  factors  (-1)  and 

(-1)  ^  of  equations  2.4.4  are  absorbed  in  the  dummy  variables 

and  £  for  the  purpose  of  this  definition  and  the  analysis  that 

follows . 


Theorem  2.4.1a 

WS  WS 

Each  of  the  component  impedances  through  Z^Q  can  be 

written  as  linear  combinations  of  the  impedance  quantity 
Z0e  (**'*;'  n'v'  S'  «)  ^ere 


1 

6  WW  f  \ 

y  ZQ  ( l ’ ,£^y,n  ,v,s,aj  dy 

0 


ZWS  ( 
Z0e  ( 


.)  = 


(2.4.6a) 


Proof 

The  proof  constitutes  finding  the  linear  combinations  of 
WS 

ZQe  ( . )  that  give  the  component  impedance  m  question 


«-(...)  =  j  ye  zf  y,  n)  dy 


ye  r  z  (jk)g 
g=0  h=0 


1  rcy'  .  i  Rcy’  1 

2  R10  2  R11  “'tYJ 


1-h 


X  z™  (..  o,  h,  g-3)dy 

1  a  1  1  1-h  1_h  rv*  1_h~h’ 

l  (jk)g  l  (i)1  h  Z  (RfX  ) 

g=0  h=0  ^  h’=0  u 


.  h*  h’ 

T 


(R^  )  U’) 


x  ZCHe+h’ )  °'h,g"3) 


(2.4. 7aa) 


1  / 
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Z2e  *  *  *  *  >  '  azt  Z0e  ("-0'  °'_1) 


««•••>  '  z  <*>»“  ''b 


h=0 


1  <Ag  ) 

h'=0  10 


,1-fc-h' 


tvi  “  1 

(ql  )  (in 


,ws 


(2.4. 7ab) 


(...e2) , 


,WS, 


«*<->  -  i  a»  zJ-9  J  V  ,1,1-h 

9=0  25  010  h=0  h'=0  2 


0  (e+h')  h,o) 
(2.4 . 7ac) ' 


(A10  > 


,  1-h-h 1 


(ACy  )  (11) 

11  T 


i  h ' 


x  (... (~2) ,  g+h ,  - 1 ) 

0 (e+h ' ) 


(2.4.7ad) 


End  pf  Proof 


Theorem  2.4.1b 

sw  sw 

Each  of  the  component  impedances  Z^f  through  Z4f  can 

be  written  as  linear  combinations  of  the  impedance  quantity 
SW 

zt\t  U'»y'*fc  ,v,s,a)  where 

1 

f  ww 

y  z”  (t' ,y' l  y,v,s,a)  dy  (2.4.6b) 

0 


The  proof  constitutes  finding  the  linear  combinations  of 
SW 

Zq^( . )  that  give  the  component  impedance  in  question.  Similar 

to  equations  2.4.7a, 


ZSW( 

Z0f( 


Proof 


Z^(...)  =■ 


1  a  1  1  i-h  1-h 

E  (3k>9  E  (i)  Z  <R^) 

J=0  h=0  z  h ' =0  lw 


1-h-h ' 


(Rjj)  Ur) 


z|*  ( . . . ) 


x  ZSW  f 

0{f+h') 


az?  ZOf  °*  °'  “1) 


0,  h,  g-3) 


(2.4. 7ba) 


(2. 4. 7bb) 


eW  1  ,  1-h  1-h  „  1-h-h'  h'  h' 

h/.0  «?»*  «n>  <V 

(2.4. 7bc) 


i«(-’  *  jo  ^  X  ^n’9’  (15’9'  io  <5,1~h 


l  (Ajj)1  n_h'  (A^)h'  U)a' 
h'=0  1  11  C 


Z0 ( f +g ' +h ' ) <•••<" 2 } - 9+h ' ~ 1 > 


( 2 . 4 . 7  bd ) 


Equation  2.4.7bd  is  dissimilar  from  2. 4. 7  ad  because  z 
is  a  function  of  n  but  not  of  y'.  This  is  easily  seen  by  writing, 
from  equations  1.4.1  and  1.4.2,  that 


Z110  zoio  Z00 0  +  aznn' 


End  of  Proof 
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Closed  form  solutions  for  are  developed  by  substituting 

equation  2.3.18  into  equation  2.4.6a. 


ZWS  f 
Z0e  (< 


)  = 


y6  Z™  ( 2. 1 ,  2  y ,  n  ,  v ,  s  ,  a )  dy 


2  (-1)0S  t  9fl  5 

a 


L 

0=  1 

2j 

sin  k2 

i 

m 

i-m 

I 

i 

t 

m=0 

n=0 

D=0 

fM  fn 

1  (_i>m  ( 

(mj  (r 

»J  1  P  l 

OeO 


(2.4.8) 


where 


C  -  (>  (*?)'’  'U.  *1101  RfP  * 


ZA2v (s+m-n) (n+a)  ^ 


(2.4.9) 


For  the  solution  of  ZgeO'  the  ar9ument  of  the  exponential 
is  first  reduced  to  within  ir/2  by  writing 
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„WS  ,  .. 

ZOeO  =  exP<"3k  r„,„  > 


010 


.cy' 


00 


*v  a  fr010Vn+P 


roio  [  tn  ] 


1 

e 

A0  1 

rno 

y 

kri 

r010 

a+n+p 


exp(-jk(rno-r010)) 


ZA2v(s+m-n)  (n+ci)  ^R])  dy‘ 


(2.4.10 


expanding  the  exponential  in  the  approximate  series  and  then 


expanding  the  binomial  (rgio  ~ rilo^  to  give 

, ,  A 


WS  tv’ 

Z0e0  =  exP<'3kr010>  A00 


a 

r010 


010 


J 


n+o 


E  a±l  (jkr010) 


il 

E 


ii-o 


m1=0 


m, 


mi 


(-1) 


1 

* 

ve 

»0  5 

frno] 

y 

3 

l*a'j 

^roio 

a+n+p+m^ 


JA2v (s+m-n) (n+a) 


( [A] , [R])  dy 


(2.4.11 
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The  factor  containing  the  integral  is 


1 

£  ] 

^(a+n+p+mi) 

1 

ye 

0 

Ao 

*0 

— ] 

s-Hn-n 

l  V 

_R_ 

a?y' 

00 

1 

ocy* 

^0  J 

C 

1 

A 

■4 

^(n+a) 


dx  dy 


*4 


^(ptm^) 


and  thus  can  be  written  as  (^Jetp^l  (M -M^here  the 

elements  A^j  of  the  matrix  [A]  are  given,  for  0  <  i,j  <  2, 

/  .  . 

(-l)0^]1  (a;)3  aJy’  i  +  j  <  2 

otherwise  (2. 


Aij  “  < 


1 

ACy' 

*00 


and  the  elements  of  the  matrix  [R]  are  given  by 


Rij  = 


*«r' 

oo 


(-i)°ia  ’  3  R-j' 


i  +  j  <  2 


otherwise 


(2.< 


and  where 


4.12) 


by 


4.13a) 


.14a) 


/ 
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ZB2NPQMS 


yMUo  yl) 


2  2  i 

z  z  x  A .  .  y3 


'■i=0  j=0 


11 


w  ,  2  2  w 

*  Uo  jio  x  Eii  10 


dx  dy 
(2.4.15) 


For  N  =  0,  ZB2NpQMS  simplifies  to 


ZB1PQMS  (fRl) 


M 


=  y 


i~o 


Oj 


X 

p 

1  J 

( 

X 

) 

2  2  +\hQ 

Z  Z  X  R.  .  y- 

i-0  j=0  ^ 

dx  dy 
(2.4.16) 


The  techniques  for  the  calculation  of  zbipqms  and 
ZB2NPQMS  are  ^etai^ed  in  Appendix  B.l  and  B.2  respectively. 
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Combining  2,4.8  and  2.4.11, 


WS  1 

Z0e  =  2j  exP(-3kr0l0} 


»5V 


a 

r010 


2 

l 

a=l 


(-l)qS 
sin  k  i 

o 


5  .  i  m  i-m 

l  a.  (jkit0)  ill 
i=0  1  m=0  n=0  d=0 


i 

m 

|(-Dm  ji-(-i) i+n+p 

m 

n 

i  P  ;  i 

roioj  n+P 

^  rr  / 


5 

r 

ii=0 


a. 

l 


(jk£*) 


,roioh 

'  i'  ' 


"i 

z 

511^=0 


fi 


1 

mi 


(-1) 


x  ZB2v(s+m-n) (n+a)e(p+m!)  <  t A 1 ' tRl > 


(2.4.17a) 


i  it  i 

The  factor  (jkrg^g)1!  has  been  rewritten  as  ( jk *-r" )  1  ( 1 
in  the  above  equation  for  two  reasons:  first,  to  make  the  expansion 

in  the  y'  direction  appear  to  be  symmetrical  with  the  expansion  in 


the  X  direction  by  writing  it  in  the  form  l  a.  ( j  k  E.  ) 1  Z  a.  (jke,  ) i  l 

X  0^1., 


X1  _1 


second  to  anticipate  a  result  to  be  derived  in  Appendix  D  whereby 


z 

Z 

Z 

Z 

m 

i-m) 

m 

n 

P 

m  1,1 

mi 

n 

P 

of 

the 

order 

r010 

n+p 

(I'L.) 

r010 

B2. 


sums  in  the  above  expression  relatively  insensitive  to  geometry 


where  a  large  variation  in  (p)  may  be  found. 


1  ;  • 
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Similarly, 


„SW  1  ...  . 

Z0f  =  2j  exP(-3kri00) 


00 


*5V 


a  _  ,  . .  os 

r100  Z. (_1) 
o=l 


s+1 


sin  ki  o 
a 


5  i  m  i-m  fml  f  i-m 

r  ai  ( iklQ) 1  E  £  Z  Ujlnll  p 
i=0  m=0  n=0  p=0 


(-1) 


m 


1-  (-1)  i+n+P 


100 


i"+p 


i1=0  *1^ 


\m,)h  ' 


m^=0 


"1 

“l 


Z32v  (s-Hn-n)  (n+ajffp+m^) 

(2.4.17b) 


where  the  elements  A^  and  of  the  matrices  [A],[R]  are  given 

for  0  =  i,j  <  2,  by 


A. 

ID 


^  (w*.)1  c,)1  *J3 


icn 

00 


1  +  D  <  2 


otherwise 


(2.4.13b) 


'ij 


_1_ 

R00 


(-1)% 


i  +  j  <  2 


otherwise 


(2.4.14b) 


•  / 


2-27 


2 . 5  Coupling  Between  Two  Surface  Segments  Revisited 
From  sections  2.1  and  2.4, 


,SS 


<l»(n)  ZWS  (n)  dn 


Ll 


4  1  f 

=  Z  Z  (-1)£  i 

5=3  f =0  ' 


: f  ZWS  ((-1)C  I'.z)  d2  . 


(2.5.1) 


The  above  equation  has  been  derived  similar  to  equation  2.4.3b. 

The  variable  2  =  ~  is  a  dummy  variable  of  integration  and 
is  unrelated  to  the  z  used  in  the  coordinate  system. 

Equation  2.4.4a  can  be  substituted  into  equation  2.5.1  to  give 


zss 


4  1 

z  z 

5=3  f=0 


(-l)f  lK 


4  1 

z  z 

t=3  e=0 


(-De  t; 


l 


Z 

i 


l 


,0 


4 

X  z 
t=l 


(e',(-i)Ti;,(-i)5 


where 

ZSS 

tef  T  5 


(2.5.2) 


(2.5.3) 


>  t 
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Theorem  2.5.1 


Each  of  the  component  impedances  z^ef  through  can  be 

written  as  linear  combinations  of  the  impedance  quantity 

zaS£  i _ , v , s, a)  where 

Oef  v  t  5 

1 

Z0ef( . >  =  j  ^  Z7e  s.«)  dz  (2.5.4) 

0 

Proof 

The  proof  constitutes  finding  the  linear  combinations  of 
SS 

ZQef;  that  give  the  component  impedance  in  question 


,SS  _ 
!lef  "  g 


1  1  i_h  l_h  ,  ,  ,1-h-h*  .  ,  ’  .  „ 

E  (jk)g  I  (i)  l  /p^y'ri)h  n  >  (R^yri) 

=0  h=0  (2‘  h'=0  (R110  h"=0  (R1C0  R101 


<y 


h" 


Z0fe+h')  (f+h")  ( •  •  .0,h,g-3) 


(2.5.5a) 


Z2ef  "  V  Z0ef 


(2.5.5b) 
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1  i  l-h  1-h  h*  .  h'  l-h-h'  l-h-h'-h"  r  .  h-  h" 


z£f=  z  (1>  1  (Anon>  {*t>  E  (Aioo‘,) 

3e£  h=0  Z  h'=0  T  h"=0 


cy'n, 


(Aioin>  <V 


x  ^e+h')  (f+h") 


(2.5.5c) 
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1 
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(A! 
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SS 

Z0(e+h') (f+g'+h") 


{...(-2) ,  g+h,  -1) 


(2.5. 5d) 


End  of  Proof 


Closed  form  solution  for  zgef  is  developed  by  substituting 
equation  2.4.17a  into  equation  2  5.4  to  give 


ZSS  ( 
Z0ef  ' 
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SS 

OefO 


(2.5.6) 
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(2.5.7) 


'  / 


■v  <***■) 
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ss 

For  the  solution  of  zoefo'  argument  of  the  exponential  is 

first  reduced  to  within  j  by  writing 
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[A000  J 

r000  1 

n+p 


j1  t 

0 

A^y’1 

00 

>5V 

fr0101 

a+n+p+i ! 

r 

r  r  ) 

Acy’n 
>•  000  J 

^r000' 

exP  j 

. 

[roio  roooj 

x  ZB2\>  (s+m-n)  (n+a)  e  (p+mj )  ^A^’  fRl)dz 


(2.5.8) 


and  expanding  the  exponential  in  the  approximate  series  and  then 


expanding  the  binomial 
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to  give 
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(2.5.9) 
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The  factor  containing  the  integral  is 

*5  (a+n+p+ij+mj) 
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000 

00 


•jsfpHn,)  x  ^(h+a) 

f  JL.V1  UJ  dx  dy  dz 


0  00 


toCY’’ 

po 


l 

RCy‘ 

f  zf 

K00 

J 

Rcy,n 

0 

K000 

%(i1-m1+m2) 


1 

e 

0 

y  1 

nCy,n 

0  1 

000 

*5  (p+mj ) 


1 

f 

[  A  )  **V  S+m-n 

R 

J 

Rcy,n 

0 

000 

K000 

^{n+o) 


dx  dy  dz 


(2.5.10) 


and  thus  can  be  written  as 


ZC2v (s+m-n) (n+a) e (p+m2 ) f (ij-mj+m2 )  ^ » ^R]) 


where  the  elements 
are  given  by 


A.  .. 
lDk 


of  the  three  dimensional  matrix  [A] 


L-' 
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Aijk  =  TsTHT  ((-D°  (^)k  A^n  ifi  +  j+k<2 


‘000 


=  0, 


otherwise 


(2.5.11) 


the  elements  of  R.  ..  of  the  three  dimensional  matrix  [R]  are 

1  J  K 

given  by 


*ijk  =  R000n  ((_1)0  "a]1  RSn  if  i  +  3+k  1  2 


=  0 

and  where 


SC2NPQMSGH  ([A^' [Rl)  " 


otherwise 


ZG  I  R  k 

Z  2  K00kZ 

lk=0  UUK 


(2.5.12) 


r1 


M 


2  2 


i  „k 


£  £  Rn,k  y32 

j=0  k=0  °3k 


*5S 


f  2  2  2  i  j  k^N 

x  111  A  x  y3  z  x 

1  ‘  i=0  j=0  k=0  > 


i  j  k^N  P 

ii"'  t 


2  2  2 


Z  Z  Z  Rijk  x  y-^z 


i=0  j=0  k=0 


i  j  k 

V  -J  7 


dx  dy  dz 


(2.5.13) 


Techniques  for  the  computation  of  zC2npqmsgH  are  treated  an 
Appendix  C.2.  Techniques  for  the  computation  of  zciPQMSGH ^R^  ' 
being  a  special  case  of  the  above  for  N  =  0,  are  treated  in 
Appendix  C.l. 
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APPENDIX  A 


CLOSED  FORM  SOLUTION  OF  SOME  SINGLE  INTEGRALS 

A . 1  The  Solution  for  Z„._- 
_ A1PQ 

Closed  form  solution  for  the  integral 
1 

ZA1PQ  ([R3}  =  /  xP  (RO4))^0  dx 

where 

R(x)  =  Rg  +  x  +  I^x2 

is  derived  here  for  various  integer  values  of  P  and  Q.  In  this 
derivation,  it  is  known  that  P  >_  0,  RQ  =  1  and  4R2  ■>  R^. 

For  Q  =  0, 

ZA1P0  (fR^)  =  P+1  (A. 1.1} 

ZA1PQ  =  ZA1P  (Q-2)  +  R1  ZAl (P+5 ( Q— 2 )  +  R2  ZA1 (P+2)  (Q-2) 

(A. 1.2) 

This  equation  is  used  as  a  recurrence  relation  for  the  cal¬ 
culation  of  for  positive  values  of  Q. 


—  J 


For  negative  values  of  Q,  a  recursion  relation  to  reduce  the 
value  of  P  is  derived  first: 

Integrating  by  parts. 


2  =  x  1  r  *SQ  I1  Q 

A1PQ  P+1  I  2  (P+1) 

0 


xP+l  R*50"1  (Rl+2R2x)  dx 


P+1  ,n  ,1  QR 


:Q  1  w"l 

Jn  "  2‘p+ 


(P+1)  ZA1(P+1) (Q-2)  ~  (P+1)  ZA1 (P+2) (Q- 2) 


(A. 1.3) 


Subtracting  A. 1.3  from  A. 1.2  and  rearranging  gives 


R2  +  P+l)  2 
(X  +  2 (P+l))  R1 


Al (P+2) (Q-2)  P+l 


-  R*50  1 

"  p  TT  R  J 


1  +  2 (P+l))  R12A1(P+1) (Q-2)  "  ZAlP(Q-2) 


(A. 1.4) 


Changing  all  P  to  P-2,  Q  to  Q+2,  for  P+Q+1^0, 


XP-1RW1  il  (2p+0)  R1 
ZA1PQ  _  (P+Q+1)R2  J0  '  (P+Q+1)2R2  ZA1(P-1)Q 


(P-1)  7 

(P+Q+1)R2  A1(P-2)Q 

If  P  =  1,  and  Q  ^  -2,  this  gives 


(A. 1.5) 


A-3 


j^q+i.  y 

ZA11Q  =  (Q+2)R2J 
If  P  >  1  and  Q?  - (P+1) 


ZA10Q 


(A. 1.6) 


ZA1PQ  (P+Q+1)R2  [(R(1|!>Q  1  ’  ZA1  (P-1)  Q  - 

(P_1)  ZA1(P-2)q]  (A. 1.7) 

For  P  +  Q  +  1  =  0,  equation  A. 1.3  rewritten  as 

ZA1PQ  =  R^  ^ZA1 (P-2) (Q+2)  “  ZAl(P-2)Q  '  R1  ZAl(P-l)Q 

(A. 1.8) 

can  be  used  recursively  until  P  =  0  or  1 .  Z,.„.  and  z,  ,  ,  , 

A10(-l)  All (-2) 

are  evaluated  separately  later  in  this  Appendix. 

Special  Case 

If  Rj  <<:  1  recursive  relationships  derived  above  accum¬ 

ulate  errors.  For  that  case,  equation  A. 1.7  rewritten  as 

ZA1PQ  -  pTT  (<*<l>>'i0+1  -  fp+iiQ+2)  ZAUPtUQ 

(P+Q+3)R2  2A1[p+2)Qj  (A. 1.9) 

along  with  the  approximation  that  for  R2  <<  1,  large  p 

ZAiPn  =  (R(D  )i{(Q)  (A. 1.10) 

Aipg  p^rj 

is  used  for  the  evaluation  of  ZAlpQ  for  all  P  and  Q. 

End  of  Special  Case 


/ 


A-4 


The  base  values  to  be  used  for  recursion  relations  A. 1.5. 
through  A. 1.8.  are  given  in  tables  of  integrals  (Gradshteyn  and 
Ryzhik  ,  1965)  and  are  reproduced  here. 


ZAlO(-3) 


(R2  +|  R1  )  R_!?(l) 

R1  +  R2 
r^Tr^+r^+TT 


if  MO 

if  A  =  0 


(A. 1.11) 


whe  re 


(A. 1.12) 


ZAlO(-2) 


tan-1 


if  A > 0 


Rl+R2 


if  A  =  0 


(A. 1.13) 


ZAll(-2) 


Jin  R  (1) 


1  \ 

2  R_ 


ZAlO(-2) 


(A. 1.14) 


•  / 
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A. 2  The  Solution  for  Z 


A2NPQ 


The  integral  zA2NPq  is  evaluated  here  given  that 
N  >-2,  P  >  0,  Q  >  -N-3  AQ=1,  Rq  =  1 »  4A2  >  A  4R2  >  R  2 

Lvt 

a  f  r  »  i  r  rs  1 1  _  f*  II  X.  I  1  ?°  1  I  P  II  X  1  1"°!^° 

A2NPQ 


f [A] , [R]1  = 

fl  ' 

'  1  ' 

T 

> 

o 

r 

p 

i  ’ 

T 

Ro] ' 

X 

Ai 

X 

X 

Ri 

(.  J 

0 

.  X2 

^  2 , 

X2  . 

r2  J 

dx 

(A. 2.1) 


is  reduced  by  using  the  recursion  relations 


ZA2NPQ  ZA2(N-2)PQ  +  A1ZA2 (M-2) (P+1) Q  +  A2  ZA2 (N-2) (P+2) Q 


(A. 2. 2) 


and 


ZA2NPQ  ~  ZA2NP(Q-2)  +  R1  ZA2N (P+1) (Q-2)  +  R2  ZA2N (P+2) (Q-2) 


(A. 2. 3) 


until  N  and  Q  are  either  zero  or  negative.  If  N  =  0, 


ZA20PQ  ZA1PQ 


(A. 2.4) 


If  Q  =  0 


ZA2NP0  ZA1PN  (fAl)- 


(A. 2. 5) 


If  both  N  and  Q  are  negative,  a  recursion  relation  on  P 
is  derived  by  integrating  A. 2.1  by  parts  to  give 


'  / 


A- 7 


ZA2NPQ 


P+1 

X 

P+1 


(A) 


*SN 


<R) 


3*0 


-il 

I 

J0 


AqR1Q  +  RqA^ 

2 (P+1)  ZA2(N-2) (P+1) (Q-2) 


AjRj^N+Q)  +  2A0R2Q  +  2RqA2N 

2 (P+1)  ZA2(N-2) (P+2) (Q-2) 

A1R2N  +  R^A2Q  +  2A2R1N  +  2R2A1Q 

2(P+1)  ZA2(N-2) (P+3) (Q-2) 


a2r2(n+q) 

P+I  ZA2(N-2) (P+4) (Q-2) 


Equations  A. 2. 2  and  A. 2. 3  have  been  used  in  deriving  the 
above  equation  wherever  necessary.  Another  equation  is  obtained 
by  applying  equation  A. 2. 3  to  each  of  the  terms  on  the  right 
side  of  equation  A. 2. 2. 


ZA2NPQ  =  A0R0  ZA2 (N-2) P (Q-2) 

+  (Aj^Rq  +  RjAq)  za2(n-2)  (P+1)  (Q-2) 

+  (A2Rq  +  +  AqR2)  za2(n-2)  (P+2)  (Q-2) 

+  (A2Rl  +  A1R2)  za2(N-2) (P+3) (Q-2) 

+  A2R2  ZA2(N-2) (P+4) (Q-2) 


'  / 


*  *:iS**M 


■  • 


A-9 


A0R 


A2R2 


;  (p*j«-2)  * 

~ N+P+Q+T 


Vo 

A2R2 


^jN+P-2) 


JA2N(P-3)Q 


(P-3)  A0Rq 


(N+P+Q+l)  A^  A2N (P-4) Q 


(A. 2. 9) 


This  recursive  relationship  is  used  while  P  >  3  and 
N+P+Q+l  ^  0.  ZA2NPQ  ^or  p  =  0»l'2  is  evaluated  later  in  this 
Appendix . 

For  N+P+Q+l  =  0,  equation  A. 2. 7  rewritten  as 


ZA2NPQ  A2R2  ZA2(N+2) (P-4) (Q+2) 


ZA2N(P-1)Q 


A1R1  .  ^0 
Aj  R  2  A2 


ZA2N(P-2)Q 


A1R0 


A2R2 


A0R1 


A2r2 


ZA2N(P-3)Q 


A°2R°2  ZA2N(P-4)Q  (A. 2. 10) 

can  be  used  recursively  until  P  is  0  or  1  or  2  or  3.  However, 
for  the  given  range  of  values,  this  equation  need  never  be 
invoked. 


/ 


A-ll 


ZA2NPQ  evaluated  for  P= 0,1,2  by  writing 


z  ,  A  ha  R  hQ 

“A2NPQ  A2  K2 


2  2 
R, 


(A. 2. 13) 


Ai  Ri  ‘'l 

If  t—  =  — ,  the  substitution  t  =  2x+=—  is  used  to  give 
A2  R2  R2 


!A2NPQ  *  2N+P+Q+1 


r 

4A2-A2 
4.2  ,  *  1 

»SN 

L.^l 

P 

*sQ 

Ri 

t  + 

A22 

i 

r2j 

"2'  J 

r2 


dt 


.W-P+O+l 


P 

(p)  (  Rif~s  7 

R1 

R1 

Z 

s=0 

lsJ  r2 J  ZA3NsQ 

R2 

'2  +  ^'  Var 

> 

(A. 2. 14) 


f  *sN  ’?*> 

(t0,t1,A1,A2)  =  (t2  +  Aj)  t  (t2+A2)  dt  (A. 2. 15 


*5Q 


where  Z 


A3NsQ 


and  where  A ,  = 


4A2  -  Aj' 


A  2  2 


ar  - 


4R2  -  Rj2 


R2' 


Ai  R^ 

If  7—^0—,  the  substitution 
A2  R2 


x  =  a  +  b  B  (t) 


where  B{t)  =  is  the  bilinear  transformation  on  t  and 


•  / 


A-12 


where 


_  A2  '  R2 

a  -  t— = - : — =—  ;  b  = 


A1R2  ~  A2R1 


,2  Rl~  A1 


AlVA2Rl 


a2+f  +  air 

A2  A2 


R2 


(A. 2. 16) 


The  above  substitution*  gives 


A-13 


Z  =  a  HQ  (2b)»sN+  JsQ  +  1 

ZA2NPQ  A2  R2  UD) 


Al] 

*$N 

Rll 

a+b  +  2A^ 

a  +  b  +  2iT 
2R2; 

(a  +  b)1 


P 

x  Z 
s=0 


:)(‘©-r 


A4N(s-N-P-0-2)Q 


ta  2A„ 


,B 


b  1 


a  + 


2Iv 


(A. 2. 17) 


where 


ZA4NPQ  (t0'tl'Al'A2) 


*sN  p  ,  >5Q 

(t2  +  Ax)  (t+lT  (t2  +  A2)  dt 


to 


(A. 2. 18) 


Beginning  of  ZA2(-1)P(-1)  evaluati°n 


A2  (-1) P (-1) 


A  ]  Ri 

is  evaluated  using  equation  (A.  2 .14)  if  —  =  -=• 

a-  r2 

A1  R1 

and  equation  (A. 2. 17)  if  —  ?  ^  ,  or  equation  (A.2.10)  if  P=1 


End  of  Za2(-1)p(-1)  evaluati°n 


_  2 


A-14 


Beginning  of  ZA2(.2)p(_1) 


evaluation 


Beginning  of  ZA2(_2)2(_1} 


evaluation 


A2 (-2)2 (-1)  "  7  _  4 

A  ^  ■K  2 


i  f  1 


1  A1  y)  1  R1 

-7—  +-T—  X+X2  =-  +  —  X  +  X2 

.  A2  A2  i  R2  r2 


dx 


2  R2**  0  [l+Jlx  +  x2] 

^2  R2 


dx 


^  Z 

A 2  A2(-2)l(-l) 


A2  ZA2 (-2) 0 (-1) 


or 

1  A1  1 
ZA2(-2)2(-l)  =  ZA10{-1)  (  fR])  "  T~2  ZA2  (-2)  1  (-1)  "  A^  ZA2(-2)0(-l) 


End  of  Z 


A2 (-2) 2 (-1) 


evaluation 


ZA2(-2)P(-1) 


for  P=0  or  1  is  evaluated  using  equation  (A. 2. 14)  if 


(A. 2. 17)  if 


,  ^1 

A2  R2  ' 


End  of  Z 


A2 (-2) P (-1) 


evaluation 


'  / 
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A. 3  The  Evaluation  Procedure  for  Z 


A3NPQ 


A3NPQ 


(t2  +  Ax)  ^  tP  (t2  +  dt 


HQ 


Beginning  of  ZA3(_^)P(_U  evaluation 
If  a2  <  ap  we  write 

ZA3(-1)P(-1)  (t0'tl'al'a2^  =  ZA3(-1)P(-1)  (fcO'tl'a2'0‘l^ 

Prom  this  point  on  it  is  assumed  that  a2  >  aj 

k  _k  _L 

Let  t  =  a i  tan©;  ©^  =  arctan  aj  ;  0Q  = arctan  aj  t^ 


a!*5*1’*  (tan©)1,  sec20  a^d0 


JA3(-1)P(-1) 


9n  (« l « 2)  **  sec2© 


—  — 

—  sm20  +  cos2© 
“2 


r*  1 


o2 


tan  0  d0 


sin29) 


/ an  a  1  ai 

Letting  A  (6)=  v  1  -  1-— Isin2©  ,  k2  =  1  -  —  ,  k'2  = 


a  2 


«2 


zA3 (-D 0 (-1)  ~  " |fC®i»M  -  *(e0,k) 


(GU  2.584.1) 
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JA3 (-1) 1 (-1) 


■“11 

**  1 

,a  2  J 

2k  ' 

In 

A  (e  i)  +k"  A(e0)-k''> 


(G  &  R  2.584.73) 


fA(6l)+k'  A[e0)-k^ 

AK)+k'  Afe^-k' 


A3  (- L) 2 ( -1)  azbk.2 


1  1  A(e1)tanei  -  E(eL,k)  -  A(9o)tan0Q 

+  E(e0,k) 


(G  &  R  2.584.90) 


-  M* 


a (e^) tane1  -  a (eQ) taneQ 
-E(e1,k)  +  E(e0,k) 


Where  F(0,k)  and  E(e,k)  are  elliptic  integrals. 
When  ai  =  a 2'  these  simplify  to 


ZA3(-1)0(-1) 


K-0O) 


ZA3(-1)1(-1) 


t  2  +  aj 
=  log - 

V  +  ai 


ZA3 (-1)2 (-1) 


^l-1^  “  ai  ZA3  (-1)0  (-1) 


End  of  Z, , .  .  .  Evaluation 

A3 (-1) P( -1) 
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Beginning  of  (—2) P(— 1)  Evaluati°n 


ZA3(-2)P(-1)  evaluate<^  using  the  substitution  t  =  ct2 
I5P  f®1  cos8  tanP6  d6 


A3 (-2) P (-1) 


=  a  2 


aj  cos2©  +  a  2  sin2© 


-w  -t, 

where  0Q  =  arctan  a2  t^;  0^  =  arctan  a2  t^. 


For  P  =  0,  the  substitution  u  =  sin©  gives 


ZA3(-2)0(-l)  =  ZA3  (sinV  Sin9l'  a2_3l' 


where 


ZA3  ^U0 '  V  b'  3)  = 


U1 


du 


UQ 


a  +  bu2 


'ab 


b,  **  1 

jarctan  ( 

!)  UlJ 

arctan 


d)V 


if  ab  > 


=  iK-  uo^ 


if  b  = 


1  U1  ~  u0 
b  uqUi 


2 (-ab) 


In 


a  +  (-ab)^u,  a  -  (-ab)^u. 


a+(-ab)^uQ  a-C-abJ^Uj 


if  a  = 


if  ab 


For  P  =  l,  the  substitution 


u  =  cos 9  gives 


'  / 


tan6 


0, 

0 

0 

<  0 
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cose 


0' 


cose.^ ,  ( 


«r°2 


End  of  Z_,,  Evaluation 

A3 (-2) P (-1) _ 


f 
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A. 4  The  Evaluation  Procedure  for  Z 


A4NPQ 


A  procedure  for  evaluating 


ZMNPQ(t0'ti'al'“2)  =  (t2  +  a  i )  (t  +  l)*'  (t2+a2)S5Udt 


HQ 


is  discussed  in  this  section  for  P  =  1,0, -1,-2.  If  P  =  0 


ZA4N0Q  ( . )  ZA3N0Q( . * 


If  P  =  1, 


ZA4N1Q( . }  ZA3N0Q( . }  +  ZA3N1Q( . * 

ZA4NPQ< . )  for  negative  values  of  P  is  invoked  only  when 

N  =  -l,  Q  =  -l.  In  that  case,  if  a2  <  aj,  we  write 

ZA4(-1)P(-1)  =  ZA4(-1)P(-1)  ^fc0 '  ^1  * 0  2  » a  l ) 

From  this  point  on  it  is  assumed  that  a2  >  aj.  Let 

H  H  H 

t  =  a i  tan6 ;  t0=O!  tanSg;  t^  =  on  tan9^. 


(l  +  o i **  tang) P  de 


a  4 (-1) p (-1)  =  H 

r  «n 

"5 

“2  9 

1-1 - sin29 

0 

l  “zj 

1  / 


A- 20 


Letting  A  (9)  =  /l  -  k2  sin20  where  k2  =  1 - -  ,  k ' 2  =  — 

«2  02 


we  write 


«i 


H  p 


A4  (-1)  P  (-1) 

a  2 


(c  +  tan9)  d9 

rp) 


where  c  =  aj 


-h 


M(-D  (-l)l-D  (l+ai)a^ 


F(9l,k)  -  F(e0,k)] 


+ 


“i 

(l  +  ai)  02*5 


IT  (9l ,  -  (1+aj  )  ,k)  -  TT  (eQ  (1+0!  )  , 


_ 1 

2/(1+<m)  ( 1+a 2 ) 


In 


d4—)*5  +  U+^)!5  A(6i) 

(l+i)>5-(l+J-)li  A  (6 1 ) 


<l+jL>  -  g+A.)  Mep) 


/ 


—  i 


Also  from  tables  of  integrals,  (G  &  R  2.591.2) 


a  (ej) 


A4(-l)  (-2)  (-1)  ojoT5  (c  +  tan  8^  cos^  (c  +  tan60)  oos2e0  <1+aiH1+<*2) 


012  i  “l  2  .*5 

(1  +  o\)  (1  +  a2)  ^5  ^+a2  +  “2  ZA4(-1)  (-1)  (-1) 


+  tanSj  A  ($i)  -  tan0Q  A (0 0 ) 


—  F(0  j  ,k)  +  —  F  (0o,k) 
a2  a.2 


-  E(0 j ,k)  +  E(0o,k) 


End  of  Z„.,  ...  ,  , ,  Evaluation 

_ A4 (-1) (-2) (-1)  _ 


The  elliptic  integrals  used  in  this  evaluation  are 
computed  using  standard  routines  from  the  SSP  library. 


V*v  *-  -*■»'**►*.*,.. 


A -22 


A.  5  The  Solution  For  ZA ( .  . . ) 

_ A,"NQMS _ 

,  1 

ZA5NQMS  (CA3'  C  B3  '  C c  3  >  =  J0  (A(x))!sN(B(x))J5Q  xM  (C (x) )hS  dx 

(A. 5.1) 

is  invoked  only  for  N=0  or  2.  If  N=0, 

ZA50QMS  =  Za2QMS  (  CB3'  tc])*  (A. 5. 2) 

If  N=  2 

2 

2A52QMS  =  AiZA2Q(M+i)S  (t33'[c])  (A. 5. 3) 


/ 
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APPENDIX  B 

CLOSED  FORM  SOLUTION  FOR  SOME  DOUBLE  INTEGRALS  ARISING  IN 
THE  COMPUTATION  OF  MUTUAL  IMPEDANCE  QUANTITIES 


Notation: 

2  2-i 

R  =  R(x,y)  =  R  (y)  =  R  (x)  =  Z  l  R.  .  x1  y 3 
x  y  i=0  j=0  13 

Writing  R(x,y) ,  as  the  trinomial  R^(x)  implies  the  form 


R  (x)  =  r 
y  i=0 


2-i 

i  Rij  y3 
j=0 


Writing  R(x,y)  as  the  trinomial  R^(y)  implies  the  form 


*x<y>  -  * 

3=0 


2-j  i 

I  Rij  x 
i=0 


[Ry] 

is  the  vector 

+  R0i1'*R02i'2){H10  +  Ruy) 

R„„ 

20 

l*V 

is  the  vector 

0 

0 

+  RjgX  +  R20x2)  (r01  +  Rnx) 

R02 

Similar  notation  is  used  with  A(x,y).  Furthermore,  [R  Q]  rep- 
rcoents  the  vector  (Rqq  Rq1  Rq2)  whereas  [ Ry=  0 ^  rePresents  the 
vector  (Rqo  R1q  R2q)  ,  etc.  The  trinomial  Rx=0(y)  is  written 
Rq (y)  for  brevity. 


'  / 
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B.l  The  Computation  of  zb1PqMS 


By  definition, 
1 


JB1PQMS 


/  yM  (R0(y))lsS 


/  xP  (R(x,y))l5Q  dx  dy 


(B.1.1) 


_  n 

Writing  (R)  =  (R)  R2  it  is  easily  proved  that 


2  2-i 

ZB1PQMS  “  LlQ  *  Rij  ZBl(P+i) (Q-2) (M+j)S 


(B.l. 2) 


and 

Z 


2 

z 

j-0 


B1PQMS  R0j  ZB1  PQ  (M+ j ) (S-2) 


(B.l. 3) 


Equation  B.l. 2  is  used  recursively  until  Q  is  zero  or 
negative.  Equation  B.l. 3  is  used  recursively  until  S  is  zero 
or  -1.  A  recursion  relation  for  reducing  the  value  of  M  and 
P  is  derived  as  follows: 

Integrating  the  inner  integral  in  equation  B.1.1  by  parts, 


B1PQMS  ~  pTT 


P+1 


yM  (Rytx))*50  dy 


-il 


x=0 


2  2-i 


P+1  . 


i=0  j=0 


i  Q  R.  .  Z 


ij  Bl(P+i)  (Q-2)(M+j)S  (B.1.4) 


'  / 


Subtracting  equation  B.1.2  from  B.1.4 


°  =  xP+1  ZA2QMS  ( tRy (x) ] , [Ry (0) ]) 


x=0 


j!0  Q  +  P+1^  Rij  ZBl(P+i)  (Q-2)  (M+j)S 


(B.1.5) 


In  equation  B.1.5,  the  sum  of  the  subscripts  in  the  P  and 
M  positions  of  ZD1  is  (P+M  +  i  +  j).  This  is  highest 

Oil* • ( 

when  i  +  j  =  2.  Separating  all  terms  that  meet  this  criterion 


I  (P  +  >5  iQ  +  1)  R..,  Z 


i=0 


i ( 2 — i )  Bl(P+i)  (Q-2)  <M+2-i)S 


-il 


xP+1  Z  t  ) 

A2QMS  1  *  * ' 


x=0 


~i=0  j=0  (P  +  ,S  lQ+L)  Rij  ZB1  (P+i)  (Q-2)  (M+ j)s 


or  changing  P  to  P-2,  Q  to  Q+2, 


E  (P  +  %iQ  +  i-l)  R. Z. 


i=0 


i ( 2 — i )  Bl(P-(2-i) )Q (M+2-i) S 


=  xP_1  Z  ( 

x  A2(Q+2)MS' 


">] 


1 

x=0 


i=0  j=0  (P  +  J5iQ  +  i_1)  Rij  ZBl(P-(2-i))Q(Mfj)S 


(B.1.6) 


_ i 
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The  implications  of  this  recurrence  relation  are  shown  in 
Figure  B.1.1.  For  a  given  value  of  Q  and  S,  different 
ZB1PQMS  va^ues  maP  onto  a  grid  in  the  P-M  plane.  The  points 
of  the  grid  with  two  circles  around  them  (A  and  G)  are  the  given 
values  of  PQ  and  Mg.  When  PQ  ?  1  (Point  A)  equation  B.1.6 
gives  a  linear  combination  of  Points  A,  B  and  C  in  terms  of  D, 

E  and  F.  When  PQ  =  1  (Point  G) ,  a  linear  combination  of 
points  G  and  H  is  given  in  terms  of  Point  J.  In  the  grid  in 
the  P-M  plane  as  shown  in  Figure  B.1.2,  it  is  assumed  that  all 
ZBl....for  which  P  +  M  <  n  -  1  are  known.  Let  Zpn  4  ZBlpQ(n_p)s  # 
where  all  Zpn  (P  =  0  to  n)  are  to  be  determined.  The  application 
of  equation  B.1.6  for  all  Zpn  (P  =  1  to  n)  gives  n  linear 
equations  to  be  solved  simultaneously.  Another  equation  is 
generated  as  follows: 


Integrating  B.1.1  by  parts. 


BlPQMS  M+l 


M+l 


**S 


y  (Ro<y>>  zAipq  <[Rx<y>]> 


1  y=0 


2  2-i 

"  ±i0  ji0  h  Qj  Rij  ZBl(P+i)  (Q-2)  (M+j)S 
2 

j=0  558:5  R°3  ZBlPQ(M+j)  ( S  —  2 ) 


(B.1.7) 


'  / 


Being  Calculated 
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Substituting  equation  B.1.2  into  the  last  term  and  equation 
B.1.3  into  the  second  term,  equation  B.1.7  becomes 


B1PQMS 


1 

M+l 


M+l 


(R0<y)> 


HS 


JA1PQ 


( [Rx(y> 1 ) 


Jy=0 


,  2  2-j  2 

"  rriT  E  2  £  (%  jQ  +*5  kS) 

”  j=0  i=0  k=0 

x  ZBl(P+i) (Q-2) (M+j+k) (S-2) 


Rij  R0k 


(B.1.8) 


Combining  equations  B.1.2  and  B.1.3,  another  equation  for 
ZB1PQMS  13  obtained* 

2  2-j  2 

ZB1PQMS  =.^Q  i^Q  k^Q  Rij  R0k  ZBl(P+i) (Q-2) (M+j+k) (S-2)  (B 


Subtracting  B.1.9  from  B.1.8  and  multiplying  throughout  by  (M+l), 

T  1 


0  =  y,M+1  (Rq{y))HS  za1pq  {  [rx  (y)  ] ) 


-ly=0 


2  2-j  2 

E  E  E  (*5  j  Q  +M+1  +  *5  kS )  R.  .  R_,  Z 
j=0  i=0  k=0  Bl(P+i) (Q-2) (M+j+k) (S-2) 


(B.1.10) 


B-8 


In  equation  B. 1.10,  the  sum  of  the  subscripts  in  the  P  and  M 
positions  ZQ1  is  (P+M+i+j+k) •  This  is  highest  when  k=2 

■71  »  •  •  • 

and  i  +  j  =  2.  Separating  all  terms  that  meet  this  criterion, 

2 

jf0  (,S  jQ  +  M+1+S)  R  ( 2- j )  j  *02  ZB1(P  +  2  -  j)  (Q  -  2)  (M  +2  +  j)  (S  -  2) 

1 

-M+1(v^sW°'  0 

-1  y=0 

2  2-j  min(3-i-j,2) 

-  E  E  t  {*5  jQ  +  M  +  1  +  >5  kS  )  R.  .  R.. 

j=0  i=0  k=0  13 

x  ZB1(P  +  i)  (Q  -  2)(M+j+k)  (S  -  2) 
or  changing  Q  to  Q  +  2,  MtoM-4,  S  to  S  +  ? , 


■  l0  (J5iQ  +M  +  S  ~  1  +  3)  R(2-j)j  R02  ZBl(P+2  -  j)Q(M~2  +  j)  S 


1 

y=0 


2  2-j  min (3-i- j , 2) 

-  E  E  E  ih  jQ  +M  +  h  ks  -3  +  j  +k)  R.  .  Rn. 

j _ 0  i-Q  K=0  ^  ^ * 


ZB1(P  +  i) Q (M  -  4  +  j  +k)  S 


(B.1.11) 


B-9 


The  implications  of  this  recursion  relation  are  shown  in 
Figure  B.1.3.  For  a  given  value  of  Q  and  S,  different 
values  map  onto  a  grid  in  the  PM  plane.  The  point  with  the 
circle  on  it  is  a  given  value  of  PQ  and  MQ.  The  left-hand 
side  of  equation  B.1.11  is  a  linear  combination  of  the  points  in  the 
shading.  The  right-hand  side  maps  the  points  in  the  dashed 
area.  Thus  equation  B.1.11  calculates  a  linear  combination 

of  the  dotted  Z  ,  quantities  in  terms  of  the  quantities  marked 
with  dashes  (which  have  M  ranging  from  Mg  -  4  to  Mg  -  1, 

P  from  Pg  to  Pg  +  2) .  It  is  also  to  be  noted  that  the  left¬ 
most  three  points  correspond  to  j=k  =  0.  Therefore,  the 
multiplier  for  those  points  in  equation  B.1.11  is  (M-3)R^gRgg. 

Thus  when  Mg  =  3,  only  the  terms  involving  M  from  0  to  2  are 
needed  for  the  computation  of  the  weighted  sum  of  the  dotted 
points.  For  the  calculation  of  Zpn,  equation  B.1.11  gives 
n-1  linear  equations,  yielding  a  total  of  2n-l  equations 
(n  linear  equations  from  B.1.6)  to  solve  for  the  n  +  1  values 
of  Zpn  The  n  +  1  equations  are  chosen  using  n  .  equations 
resulting  from  letting  P=n  down  to  1 (M  =  n  -  P)  in  equation 
B.1.6  and  1  equation  resulting  from  letting  M  =  n  (P=  0)  in 
equation  B.1.11.  This  set  has  been  chosen  so  as  to  maximize 
the  use  of  equation  B.1.6  which  is  simpler  in  form. 


I 

1 


Pictorial  Representation  of  Equation  B  .1 .11 
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The  application  of  equation  B.1.11  for  P  =  0,  M  =  n,  gives 


(n  +  S-1)  R2Q  RQ2  Z2n  +  (n +S  -  1  +  5s  (Q  +  2) )  R, ,  Rrt0  Z 


11  02  "In 


+  (n  +  S-  1 +  ( Q  +  2 ) )  R02  RQ2  ZQn  =  Sln  (B. 1.12a) 


where 

S 


In 


-  y"’3  (R0(y))'i!s+2)  zA10(Q+2)  (I*x(y))) 


2  2-1 


y=0 


Z  £ 
j=0  i=0 


3  min(3-i-j,2) 


k=0 


(h  jQ  +  n  +  %  kS  -  3  +  j  +  k)  Ri ^  RQk 


Bl(0+i)Q(n-4  +  j  +  k)  S  (B. 1.13a) 

The  application  of  equation  B.1.6  for  P  =  1  and  2  (M  =  n  -  1  and 
n  -  2  respectively)  gives 


*  (Q  +  2)  R11  Z0n  +  <Q+2’  R20  Zln  =  S2n 


(B. 1.12b) 


R02  Z0n  +  (1  +  *  <Q  +  2) )  Ru  Zln  ♦  (1  +  (Q  +  2) )  R2Q  Z2n  -  (B.  1.12c) 

where 

S2n  =  ZA2  (Q+2)  (n-1)  S  ^Ry(1)  1  '  [RY(0>  l)  "  ZAl(n-l)  (Q+S+2)  ^lRy(0) 

1  1-i 

~  i=0  j=0  *i{Q  +  2)  Rij  ZBl(l-(  2-i))Q(n-l+j)S 


(B. 1.13b) 


B-12 


and  where 

S3n  ZA2(Q+2)  (n-2)S  '  fRy(0)l) 


1  1-i 

l  Z  (1  +  \  i  (Q+2)  )  R . ..  Z 
i=0  j=0 


ij  B1  i Q (n+ j-2) S 


(B. 1.13c) 


Equations  B.1.12  are  solved  simultaneously  to  give 


z~  Z-  and  Z-  .  It  is  seen  that  the  three  equations  are  singular 


On,  “In  -  “2n- 

if  Q  =  -2  in  which  case  Z^^  and  Z^n  are  computed  numerically  and 
Z2n  is  computed  using  equation  B. 1.12a.  Once  ZQn#  Zln  and  Z2n 
are  known,  all  other  Zpn  are  computed  by  writing  equation  B.1.6 
in  the  form 


JPn  (P+Q  +  1)R 


20 


ZA2(Q+2)MS  nRy(i)]  '  [Ry(0)]) 


(P  +  *5  Q)  Ru  Z(p_1)n  -  (P-D  r02  z(p-2)n 


1-i 


-  z  Z  (P  +  H  (Q+2)  -  1)R.  .  z, 


i=0  j=0 


ij  B1 (P- ( 2-i) ) Q (n-P+  j )  S 


(B.1.14) 


Figure  B.1.4  summarizes  the  calculation  strategy. 

A  recursion  relation  that  is  simpler  than  equation  B.1.11 
can  be  derived  when  S  =  0.  Under  this  condition  the  last  term  in 
equation  B.1.7  vanishes  and,  subtracting  B.1.2  from  B.1.7, 


*  / 


Figure  B.1.4:  Strategy  for  the  Calculation  of  Z^p^s  for  Various  Values  of 
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1 

-  yM+1  zAiPo([E*<i,)l! 

-I  v  =  f 


2  2-i 

-  I  E  (>5  jQ  +M  +  1)  R.  .  Z 


i=0  j=0 


ij  Bl(P+i) (Q-2) (M+j)0  (B.l. 


Separating  all  terms  for  which  i  +  j  =  2, 


Z  (k  jQ  +M  +  1)R 


( 2- j ) j  Bl(P+2-j)  (Q-2)  (M+j)0 


=  y  ZA1PQ 


1  1-i 

E  E  (*s  jQ+  M+  1)  R.  .  Z. 


i=0  j=0 


ij  Bl(P+i) (Q-2) (M+j)0 


or,  replacing  M  by  M-2,  Q  by  Q+2  throughout. 


jio  (Js3Q  +  M-l+j)  R(2-j)j  ZBl(P+2-j)Q(M-(2-j))0 


y  ZAlP{Q+2)  ^Rx(y)  ^ 


Jy=o 


1  1-i 

j!0  (!s  3  (Q+2)  +M-1)  j  ZBl(p+i)Q(M-(  2- j ) )  0 


(B.l. 16) 
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Using  B.1.6  and  B.1.16,  all  points  on  the  P-M  grid  except 
P=M=0  become  candidates  for  recursion.  When  n=l ,  M=n ,  P=0, 
equation  B.1.16  gives 

%(Q+2)  RX1  2^+  (Q+2)  Rq2  ZQ1  =  zAlP(Q+2)  (tRx(1)l) 


ZAlP(Q+2)  (fRx(°)  1)  “  *5(0+2)  R01  Zb10Q00 


(B.1.17) 


which  is  solved  simultaneously  with  equation  8.1.12b  to  qive  ZQ 
and  Z^q.  When  n=2,  equation  B.1.16  for  M=2,  P=0  becomes 

R20  Z20  +  (1  +  *5(0+2))  RX1  Zu  +  (l+(Q+2) )  RQ2  ZQ2 

=  ZA10(Q+2)  ^  tRx(1) 1) 

1  1-i 

-  S  l  «*a  j  CQ+2)  +  1)  R..  Z  (B.1.18) 

i=0  j=0  13  BllQ:)0 

Figure  B.1.5  shows  the  P-M  grid  and  outlines  the  stragety 
for  calculating  all  zBiPQMo' 

The  numerical  calculations  are  performed  using  a  $  noint 
Gaussian  integration  (Abramowitz  and  Stegan,  formulae 
25.4.  30)  whereby 
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dy 


6 

I 

i=l 


wif(2Xi  +  2* 


where  w^  and  are  given  by  (Abramowitz  and  Stegun,  Table  25.4) 
for  various  values  of  n .  Thus, 


ZB1PQMS 


i yt> "  <VV 


S/2 


ZA1PQ 


where  y^ 


1 

2 


(x.+l) . 


(B.1.19) 


'  / 
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Applying  equation  B.2.2  to  the  second  term  and  B.2.4  to  the 
last  term  in  equation  B.2.5  and  then  subtracting  the  combination 
of  equations  B.2.2  and  B.2.4  from  it 

0  =  xP+1  ZA5NQMS  (CAx3 '[Rxl ' 

_.x=0 

2  2-i  2  2-k 

-  I  l  Z  Z  (h  iN  +P+1  +h  kQ)  A.,  R 

i=0  j  =0  k=0  2=0  13  KX' 

x  ZB2(N-2) (P+i+k) (Q-2) (M+j+*)S 

The  sum  of  the  subscripts  in  the  P  and  M  positions  is 
(i  +  j+  k  +  2+  P  +  M).  This  is  highest  when  i  +  j  =  2;  k  +  2  =  2 . 
Separating  all  terms  that  meet  this  criterion, 

2  2 

i  =0  k=0  ( **  1  ”  +  P  +  1  +  **  R  Q)  Ai  ( ■ 2~i )  (2-k)  ^2 (N-2)  (P+i-Ht)  (Q-2)  ((MMHi-*) )  f 

=  xP+1  ZA5NQMS  ( *  *  * }  ]  x=Q 

2  2-i  min{3-i-j,2)  min(3-i-j-k,2-k) 

-  Z  Z  Z  Z  (%iN  +  P+1  +  J*Q)  A.  .  R 

i=0  j=0  k=0  1=0  13  * 

X  ZB2(N-2)  (P+i+k) (Q-2) (M+j  +  2)  S 

or,  replacing  N  by  N+2,  P  by  P-4,  Q  by  Q+2,  throughout. 


-  / 
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2  2 

kf0  (  Jsi(N+2)+P“3+islt(0f2))  Ai(2-i)Rk(2-k)  ZB2N(E*- (4-i-k)  )Q(M+4-i-k)S 


=  x 


P-3 


A5(N+2)  ;Q+2)MS 


(...) 


] 


1 

x=0 


2  2-i  min(3-i-j,2)  min(3-i-j-k,2-k) 

-  Z  Z  Z  Z  (*si  (N+2)  +P-3+yc (0+2) )  A..R 

i=0  j=0  k=0  £=0  11 


x  ZB2N (P- (4-i-k) ) Q (M+j+i) S 
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The  recursion  implications  of  equation  B.2.6  are  shown  in 
Figure  B.2.1.  In  the  P-M  grid  as  shown  in  Figure  B.2.2,  it 
is  assumed  that  all  Z^2  for  which  P+M  <  n-1  are  known.  Let 

ZB2NPQ(n-P)S  where  a11  ZPn  are  to  be  determined.  The 
application  of  equation  B.2.6  for  all  Zpn  (P = 3  to  n)  gives 
(n  -  2)  equations  to  be  solved  simultaneously.  If,  as  is 

sometimes  the  case,  Z(p_1)n*  z(p-2)n'  Z(P-3)n  and  Z(P-4)n 

are  known,  equation  B.2.6  written  in  the  form 


(N  +  P  +  Q  +  1) 


2  2-i 

-  z  z 


i=0  j=Q 


1 

#20  *20  WqmS  =  XP-3  2a5(n+2)(q.2)ms  (...1 

Jx=0 

min(3-i,2)  2-k 

Z  Z  (*si  (N+2)  +P-3+*5k (Q*-2) )  A.  .  R 

k=0  £=0  13 


X  ZB2N(P- (4-i-k) )Q(Mtj+£)  S 

is  useful  for  finding  «B2NpQMS> 


(b.2.1  ) 


Being  Calculated 


If  Q  =  0,  following  a  procedure  similar  to  that  used  in  the  devel¬ 
opment  of  equation  B.1.6,  an  equation  simpler  than  B.2.6  is 
derived: 


2  i 

i=0  (2(N+2)l+P-1)Ai(2-i)  ZB2N(P-(2-i))0(M+2-i)S 

=  X  ZA2  (N+2 )  MS  (fAx^'  ]x=0 

jj0  (2(N+2)l+P~1)Aij  ZB2N(P-(2-i))0(M+j)S 

(B.2.8) 

Other  equations  are  generated  depending  on  whether  Q  and 
S  are  zero  by  integrating  the  outer  integral  in  equation 
B.2.1  by  parts 
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zb2npqms  m+i  y  tRo(y))  za2npq  (f-Ay J '  [Ry])]y 
2  2-i 

i=0  j=0  i5Nj  Aj-3  ZB2(N-2)  (P+i)Q(M+j)S 

2  2-i 

-i-n  !sQj  Rii  ZB2N(P+i)  (Q-2)  (M+j)  S 


-n  *  Sj  R°3  ZB2NPQ(M+j) (S-2) 


If  both  Q  and  S  are  zero,  equation  B.2.9  reduces  to 


7  -  1  vM+1 

B2NP0M0  “  M+l  y 


ZA2NP0  ( 


•  o]1 

Jy=0 


L  2  2-i 

M+T  iZQ  ^  Aij  ZB2  (N-2)  (P+i)  0(M+j)  0 
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Comparing  with  equation  B.2.4, 

1 

y=0 


0  =  y 


M+l 


ZA2NP0(‘ 


2  2-j 

E  E 
j=0  i=0 


(M  +  l+ijNj)  Aij  zB2(n-2)  (P+i)  0(M+j)  0 


or  separating  the  terms  for  which  i  +  j  =  2  while  replacing  N 
by  N  +  2 ,  MbyM-2  throughout 

j=0  (M"  1  +Ss(N+2)  j)  A(2-j)i2B2N(P+2-i)0(M-(2-j))0 

1 

y=0 


=  y 


M-l 


a2(N+2)P0 


(..) 


1  i-j 
E  E 
j=0  i=0 


(M  -  1  +  *5 (N+2)  j )  zB2N(P+i)0  (M-(2-j)  )0 


(B.2.10) 


Figure  B.2.3  shows  the  recursion  formula  pictorially. 

For  S  =  0,  Q/0,  equation  B.2.9  becomes,  by  application  of 
equation  B.2.2  to  the  second  term  and  equation  B.2.4  to  the 
third  term, 
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Figure 
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B2NPQM  0 
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1  M+l 

m+t  y 


2  2-i  2 

e  e  e 
i=0  j=0  k=0 


2-k 

E 
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ZA2NPQ(tAy3  JRy])]y=() 

(JjNj  +  »sQi)  A£  j  Rj^ 


ZB2(N-2)  (P+i+k)  ( Q—  2 )  (M+j+i) 0  (B.2.11) 

Combining  equations  B.2.2  and  B.2.4  and  subtracting  the 
result  from  the  above, 

0  =  Y  ZA2NPQ  { ‘  }  Jy=0 

2  2-i  2  2-k 

-  E  E  E  E  (*SN  j  +  *sQ*.  +  M  +  l)  A^  R^ 

i=0  j=0  k=0  11=0  J 

ZB2 (N-2) (P+i+k) (Q-2) (M+j+Z)  0 

or 


1  / 


0  =  y 


M+l 


ZA2NPQ ( *  * J 


Jy=0 


2  2-j  2  2-8, 


-  I  I 


j=0  i=0  8=0  k=0 


l  Z  (*5  N  j  +  JsQ8  +  M  +  l)  A.  .  R. 

—  A  1- _ A  A  J  ^  ^ 


JB2(N-2) (P+i+k) (Q-2) (M+j+8)  0 


(B.2.12) 


Separating  all  terms  for  which  i  +  j+  k  +  8  =  4. 


2  2 


j=0  z=0  ^  Nj  +  **  Ql  +  M+1^  A(2-j) jR(2-8)8*B2(N-2)  (P+4-j-8)  (Q-2)  (M+j+8)  0 


=  y 


M+l 


ZA  2NPQ 
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y=0 


2  2-j  min(3-i-j,2)  min (3-i- j , 2) -8 

l  Z  Z  Z  (*5Nj  +  HQl  +  M+l)  A.  .  R, 

j=0  i=0  8=0  k=0  +  J 


x  Z 


B2 (N-2) (P+i+k) (Q-2) (M+j+8)0 


or  replacing  N  by 


N+2 


Q  by  Q+2 ,  M  by  M-4  throughout. 


*mQ  t^Q  ('*(M  +2)j+!5<Q  +2)t+M'3)  A(2-j)  j  R(2-l)t  ZB2N(P+4-j-t)Q<M-(4-j-O)0 


=  y 


M-3 


ZA2(N+2)P(Q+2) (,,) 


y=0 


2  2-j  min(3-i-j,2)  min(3-i-j,2)-i 

ll  l  l  (*5(N+2)  j +*s(CH-2)i+M-3)  A.  .  R 

j=0  i=0  e=o  k=0  KX- 


x  ZB2N(P+i+k)Q(M-(4-j-4) )0  (B.2.13) 

Figure  B.2.4  shows  this  relation  pictorial ly. 

For  Q  =  0,  S^O,  combining  equations  B.2.3  and  B.2.4  and 
subtracting  the  result  from  equation  B.2.9  after  application  of 
equation  B.2.3  to  the  second  term  and  equation  B.2.4  to  the 
last  term, 

o  -  y"*1  (V',)5',S  W  ’•>  ]y„0 

2  2-j  2 

-  I  l  t  A.  .  Rnn,  (%  Nj  +  M+l  +  h  Sm) 
j=0  i=0  m=0  13  0m 

ZB2(N-2) (P+i) 0(M+j+m) (S-2) 


(B.2.14) 


Figure 
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Separating  all  terms  for  which  i  +  j+m  =  4. 


j=0  (ls  N  j+  M+  1  +  S)  A  (2- j)  j  R02  ZB2(N-2)  (P+2-j)  0(M+2  +  j)  (S-2) 


=  yM+1  (Rn<y)) '  za1Np  <•  > 


*5S 


A1NP 


1  1 


Jy=0 


2  2- j  min ( 3— i—  j , 2) 

-  £  Z  E  (35Nj  +  M+  l  +  )5  Sm)  A .  .  R. 

j=0  i=0  m=0  ^  0m 


x  ZB2(N-2)  {P+i)0  (M+j+m)  (S-2) 


or  replacing  N  by  N  +  2,  M  by  H-4,  SbyS  +  2  throughout. 


(Js(N+2)  j  +M-  1  +  S)  A  (2- j )  j  R02  ZB2N(P+2- j)0  (M-(2-j)  )S 
M  ,  h(S+2)  l1 

y  (R0<y))  zai(n+2)p 


-<  y  =  0 

2  2- j  min(3-i-j,2) 

-Z  Z  Z  (*5(N+2)  j  +  M-3  +  H  (S+2)  m  )  A  .  .  R 

j=0  i=0  m=0  J 


X  ZB2N(P+i)0 (M- (4-j-m) ) S 
Figure  B.2.5  shows  this  relation  pictorially. 
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figure 
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If  both  Q  and  S  are  non- zero,  subtracting  the  combination 
of  equations  B.2.2,  B.2.3  and  B.2.4  from  equation  B.2.9  after 
applying  equations  B.2.2  and  B.2.3  on  the  second,  equations  B.2.3 
and  B.2.4  on  the  third  and  equations  B.2.4  and  B.2.2  on 
the  last  terms, 

jjg  l1 

o  -  yMtl  (n0<y»  «AMP0(C\]  .  C«y]> 

Jy=0 

2  2-j  2  2-4  2  r  -| 

-l  Z  l  Z  l  MiNj  +  %Q4  +  M+l  +  >$Sm  A.  .  R.  Rn 
j=0  i=0  4=0  k=0  m=0 «-  -1  13 

x  ZB2(N-2) (P+i+k) (Q-2) (M+j+4+m) (S-2)  (B.2.16) 

Replacing  N  by  N  +  2,  Q  by  Q  +  2,  Mby  M-6 ,  S  by  S  +  2 
throughout  while  separating  the  terms  for  which  i  +  j+  k  +  4+  m  =  6, 
equation  B.2.16  becomes , 

2  2  r 

I  1  h (N+2) 
j=0  4=0  L 

X  ZB2H(P+4-i-j-4)Q(M-(4-j-4) )S 

M_5  *S  (S  +  2)  -]1 

-  y  (R0<y))  za2(n+2)p(q+2)  (,-)  J 

y=0 

2  2-j  2  2-4  min(5-i-j-k-4,2)  r  j 

_  Z  Z  Z  Z  Z  hs(N+2)  j  +  *s(Q+2)  4+M“5+*5(S+2)m 

j=0  i=0  4=0  k=0  m=0  *-  “* 

X  Aij  \i  R0m  ZB2N (P+i+k) Q (M-  (6- j-4-m) ) S  (B.2.17) 


j  +  *5  (Q+2)  4  +  M  -  3  +  S 


(2- j ) j  (2-4)4  02 


'  / 
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Figure  B.2.6  shows  this  relation. 

The  case  where  Q  =  S  =  0  occurs  whenever  the  far-field 
approximation  is  being  considered.  Calculation  strategies  are 
developed  here  for  the  Q  =  S  =  0  and  the  general  case  Q/0, 

S  /  0.  The  other  two  cases  Q  =  0,  S  f  0  and  Q  ^  0,  S  =  0  occur 
infrequently  enough  to  justify  not  programming  them  as  special 
cases  but  rather  to  use  the  general,  even  when  they  occur. 

Figure  B.2.7  summarizes  the  calculation  strategy  for  calculating 

ZB2NPQMS-  The  relation 

ZB2NPQMS  =  2  if1  wi  (yiJ  (Vyi>>  ZA2nPQ  ( ^Ay=y^  ' 

(B.2.18) 

is  used  for  numerical  integration.  Here,  y^  =  ^-(x^+1)  and  w^ 
and  are  given  by  Abramowitz  and  Stegun  (Table  25.4,  n=6) . 
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APPENDIX  C 


CLOSED  FORM  SOLUTION  FOR  SOME  TRIPLE  INTEGRALS  ARISING  IN 
THE  COMPUTATION  OF  MUTUAL  IMPEDANCE  QUANTITIES 

Notation: 


R  =  R(x,y,z)  =  Rx(y,z)  =  Ry(z,x)  =  Rz  (x,y)  =  Rxy(z) 


2  2-i  2-i-j  ^ 

=  R  (x)  =  R  (y)  =  Z  Z  Z  R.  x  yJZ  . 

zx  1=0  j-0  k=0 


yz 


Writing  Rx(y,z)  implies  the  form 


2  2-j 

R  (yfz)  =  Z 

j=0  k= 


-j  [2-j-k  i 

Z  Z  Rijk  x1 

=0  L  i=0  J  J 


yjzk, 


i.e.,  the  variable  appearing  in  the  subscript  is  the  most  suppressed. 
This  definition  is  consistent  with  the  statement  that  £RXJ  is  tfie 
matrix 
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000 

+  R100  X  +  R200X 

R001  +  R101X 

R002 

010 

+  R110  x 

R011 

0 

020 

0 

0 

Similar  notation  applies  to  the  subscripts  y  and  z.  Consistent 


with  this  notation 


C-2 


R0(y»z)  =  r 


2-j 


j=0  k=0 


E-  Rojk  y  z 


i  * 


and 


m 


represents  the  matrix 


R 


000 


010 


020 


001 


R 


Oil 


002 


Writing  R  as  R^z  (x)  implies  the  form 


R  (x)  = 
yz 


2  r  2- 

z  z 

i-o  |_j= 


2-i  2-1-3  v1 

z  Riik  y3z 

k=0  XJK 


and  is  consistent  with  the  statement  that 


[Ryz]  1 


is  the  vector 


r  2 
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2"j  i  k 

Z  y  2 

j=0  k=0  u:,K 


1  1-3  i  k 
E  E  Rlik  YDzk 
j=0  k=0 


'200 


Consistent  with  the  above  definitions,  |R 


[ *'  (y=0)  (z 


-3 


implies  the 


R000 

R100 

R200 


vector 
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C.l  The  Computation  Of  zciPQMSGH 


By  definition, 
ZC1PQMSGH 


^  zG  (R00(z))H/2  j  yM(RQ(y ,z)) 


S/2 


/  P  Q/2 

x  (R(x,y,z))  dxdydz 


(C.1.1) 


Writing  (.)n  =  (.)n  2  (.)2  it  is  easily  proved  that 

2  2-i  2-i-j 

ZC1PQMSGH  =  kfQ  Rijk  ZCl(p+i) (Q-2) (M+j)S(G+k) H, 

2  2-j  '  K-1*2* 


ZC1PQMSGH  klQ  R0jk  ZClPQ(M+j)  (S-2)  (G+k)H 


(C.l. 3) 


and 


ZC1PQMSGH  =  kfQ  R0Gk  ZClPQMS(G+k)(H-2) 


(C.l. 4) 


Equations  C.l. 2  through  C.1.4  are  used  as  recursion  relations  for 
desired  values  of  Q,  S  and  H  in  terms  of  (smaller)  given  values  of 
Q,  S  and  H.  Relations  for  giving  ZC1  for  desired  values  of  P,  M 
and  G  in  terms  of  given  (smaller)  values  of  P,  M  and  G  are  derived 
as  follows: 


! 


/ 


—  i 
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Integrating  the  integral  on  x  by  parts 

ZC1PQMSGH  =  ■  — —  xP+1  J  Z<3  (R00(Z))H/2  /  yM 

0  0 
(R*)(y/z))S//2  (\.(y,z))Q/2  1 

x*0 

1  2  2-i  2-i-j  . 

*P+1  Jo  ji0  kJ0  2  lQRijk  ZCl(P+i)  (Q-2)  (M+j)S(G+k)H 

(C.1.5) 

Comparing  with  equation  C.1.2, 

0  *  x  ZB2QMSGH  (fRx^'  CRx=oJJ"} 

J  x=0 

2  2-i  2-i-j  ^ 

if0  k^0  (2  iQ  +P+1)Rijk  ZCl(P+i) (Q-2) (M+j)S(G+k)H 

(C.1.6) 

Separating  the  terms  for  which  i+j+k  =  2,  which  changing  P  to  P-2, 

Q  to  Q+2 
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i=0 


j=0  (2X (Q+2)+P_1)  Rij(2-i-j)  ZC1 (P+i-2)Q(M+j)S(G42-i-j)H 
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B2 (Q+2 ) MSGH  (..) 

x=u 

1 

1-i 

1-i-j 
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i=0 

I 

j  =  0 

L 

k=0 

(  i  i(Q+2)+P-l)  Rijk 

X  ZCl(P+i-2)Q(M+j)S(G+k)H  (C.1.7) 


The  implications  of  this  recurrence  relation  are  shown  in 
Figure  C.1.1.  For  a  given  value  of  Q,S  and  H,  different  zcipqMSgh 
values  map  onto  a  grid  in  the  P-M-G  coordinate  system.  The  points 
of  the  grid  with  two  circles  around  them  (A  and  L)  are  the  given 
values  of  PQ,M0  and  GQ .  When  PQ/1,  (Point  A)  equation  C.1.7  gives  a 
linear  combination  of  points  A  through  F  in  terms  of  points  G,  H,  J  and  K. 
When  Pq=1  (Point  L) ,  a  linear  combination  of  points  L,  M  and  N  is  given 
in  terms  of  point  P.  A  plane  of  constant  P+M+G  is  shown  in  Figure  C.1.2 
Recursion  relation  C.1.7  can  be  applied  to  all  points  for  which 
P%1.  n+1  more  equations  must  be  generated  (for  P=0)  in  order  that 

all  points  in  this  plane  may  be  solved.  This  is  done  by  integrating 
the  integral  on  y  by  parts  to  give 
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I 
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Figure  C.1.1:  Implications  of  the  Recursion 
Relation  C.1.7 
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C1PQMSGH 
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2  Sj  R°jk  ZclpQ(M+j> (S-2) (G+k)H 


(C.1.8) 


Substituting  equation  C.1.3  into  the  second  term  and  equation  C.1.2 
into  the  third  term,  comparing  with  a  combination  of  equations 
C.1.2  and  C.1.3,  replacing  Q  by  Q+2,  M  by  M-4  and  S  by  S+2  through¬ 
out  and  writing 

ZB3P0HGS  <W'  [“]»  ~J  (s(y))H/2  yG  (B0(y))s/2 

0 

X  jV  (B(x,y))Q/2  dxdy  (C.1.9) 

where  [b]  is  a  matrix  whose  first  row  is  BQ  and  [a]  is  a  vector,  the 
following  equation  is  obtained 


°  Y  ZB3P  (Q+2)HG(S+2)  ( [R0C 3  '  [Ryp] 
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x  RijkR0im  ZC1 (P+i)Q(M+j+i-4)S(G+k+m)H 


*  Techniques  for  solving  ^pqhgs  •  M  >  are  similar  to  those  presented 
in  Appendix  B.2  for 
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or,  separating  the  terms  for  which  i+j+k=2,  £+m=2. 


j  V  I  ,i  M-l-§S.S-»)  ''(Z-j-kljk  »0<2-m> 
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m 


x  Z 
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(C.1.10) 

It  is  to  be  noted  that  equation  C.1.7  helps  to  write  zCiPqMSGH 


in  terms  of  other  ZC1  with  smaller  values  of  P  but  larger  .  v 

values  of  M  and  G.  Equation  C.1.10  gives  Z  .  in  terms  of 

or ..... . 

other  Zc^  with  smaller  values  of  M  but  larger  values  of  P 

and  G.  Another  recursion  relation  that  gives  Z  .  in  terms  of 

other  ZC1  with  smaller  values  of  G  but  larger  values  of  P  and 

M  is  derived  by  integrating  the  integral  on  z  in  equation  C.1.1  by 
parts  to  give 
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G+l  {z  (R00(z))H/2  ZB1PQMS  {IW] 

■J  z=Q 
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kfQ  2  Qk  Rijk  ZC1 (P+i) (Q-2) (M+j)S(G+k)H 

2"j  1 
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Substituting  equations  C.1.3  and  C.1.4  into  the  second  term, 
equations  C.1.4  and  C.1.2  into  the  third  term,  equations  C.1.2  and 
C.1.3  into  the  fourth  term,  comparing  the  result  with  a  combination 
of  equations  C.1.2,  C.1.3  and  C.1.4,  replacing  Q  by  Q+2,  S  by  S+2, 
G  by  G-6,  H  by  H+2  throughout,  the  following  equation  is  obtained: 

0=2  (R00(z))  (H+2,/2  ZBlP(Q+2)M(S+2)  ([Rz])1 
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2  2-i  2-i-j  2  2-1  2 
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i=0  j  =0  k=0  £=0  m=0  n=0 

(|(Q+2)k+|(S+2)m+G-5+|(H+2)n) 

Rijk  R0£mR00n  ZC1 (P+i) Q (M+ j+£ ) S (G+k+m+n-6) H  (C.1.12) 

or,  separating  the  terms  for  which  i+j+k=2,  £+m~2,  n=2 
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(C. 1 .13) 

As  suggested  in  the  discussion  on  Figure  C.1.2,  equations  C.1.7, 
C.1.10  and  C.1.13  can  be  used  for  various  values  of  P,  M  and  G  in 
order  to  generate  enough  equations  to  solve  for  all  points  on  the 
constant  P+M+G  plane  simultaneously.  Figure  C.1.3  shows  the  left- 
hand  side  points  of  equation  C.1.7  (i.e.r  each  application  of 
equation  C.1.7  yields  a  linear  combination  of  a  cluster  such  as  is 
shown  in  Figure  C.1.3).  Similarly  Figure  C.1.4  shows  the  left-hand 
side  points  of  equation  C.1.10  and  Figure  C.1.5  shows  the  left-hand 
side  points  of  equation  C.1.13.  At  least  one  of  the  equations  is 
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Figure  C.1.4:  Recursion  Implications  of 

Equation  C.1.10 
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Figure  C.1.5:  Recursion  Implications  of 

Equation  C.1.13 
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applicable  at  each  of  the  points  and  is  solved  simultaneously  to 
give  the  value  of  Zc^  at  each  of  the  points. 


The  case  where  Q=S=H=0  is  of  special  importance  because  it  is 
involved  for  the  far  field  case  (see  Appendix  D  for  details) . 

The  three  integrals  in  the  definition  of  z C1  decouple  in  this 

case  and 


(C.1.14) 
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The  Computation  of  ZC2NPQMSGH 


By  definition, 
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(C.2. 1) 
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it  is  easily  proved  that 
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0 jk  ZC2NPQ(M+j) (S-2) (G+k) H 
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R00k  ZC2NPQMS (G+k) (H-2)  (C.2.5) 


Equations  C.2.-2  through  C.2.4  are  used  as  recursion  relations 
for  desired  values  of  N,  Q,  S  and  H  in  terms  of  (smaller)  given 
values  of  N,  Q,  S  and  H.  Relations  for  giving  f°r 


desired  values  of  P,  M  and  G  are  derived  as  follows: 
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Integrating 

equation  C.2 

. 1  by  parts 

dummy  variable  x 

and  applying 

equations  i 
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(C. 2.6) 


Comparing  with  a  combination  of  equations  C.2.2  and  C.2.3  and  chang¬ 
ing  N  to  N+2,  P  to  P-4,  Q  to  Q+2  throughout. 


P-3 

0  =  X  ZB3(N+2) (Q+2) NSGH  1 


M-  M-M>L 
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(C.2.8) * 


Separating  ali  terms  in  equation  C.2.~  for  which  i+j+k=2  and 
£+m+n=2,  the  following  recursive  relation  is  obtained: 
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(C.2.9) 


*  Techniques  for  the  solution  of  Zb3nqmsgh([a] ,[b], [C])  are 

similar  to  those  presented  in  Appendix  B.2  for  ZB2NpQMS  (fA] ,  j 
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If  equation  C.2.9  is  written  with  N**Nq ,  P=PQ,  etc., 
each  of  the  terms  on  the  left-hand  side  has  P+M+G=Pq+Mq+Gq. 

The  terms  on  the  right-hand  side  have  P+M+G=P0+M0+GQ-n  where 
n  is  1,  2,  3,  or  4.  This  statement  implies  that  none  of  the 
points  on  the  right-hand  side  of  equation  C.2.9  fall  on  the 
(Pq+Mq+Gq=  constant)  plane.  Figure  C.2.1  shows  a  plane  of 
constant  Pq+Mq+Gq.  Recursion  relation  C.2.9  may  be  applied 
to  all  points  on  this  plane  for  which  Pq»3.  Equations  that 
may  be  applied  for  smaller  values  of  PQ  are  developed  next. 
Integrating  equation  C.2.1  by  parts  with  respect  to  the  dummy 
variable  y  and  applying  two  of  the  equations  C.2.2  through  C.2.4 
to  each  term  except  the  first, 
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Comparing  equation  C.2.10  with  a  combination  of  equations 
C.2.2  through  C.2.4  ,  changing  N  to  (N+2) ,  Q  to  (Q+2) ,  M  to 
(M-6) ,  S  to  (S+2)  throughout,  and  separating  all  terms  for  which 
i+j+k=2,  £+m+n=2,  and  p+q=2,  the  following  recursion  relation  is 


obtained: 
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(C.2.12)* 


*  Techniques  for  the  solution  of  zB4NP^iGS  are  similar  to  those  presented 


in  Appendix  B.2  for  . 
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Here,  [  A]  and  [B]  are  matrices,  Bq  is  the  first  row 
of  the  matrix  [s]  and  [c]  is  a  vector. 

Integrating  equation  C.2.1  by  parts  with  respect  to  the 
dummy  variable  z  and  applying  three  of  equations  C.2.2  through 
C.2.5  to  each  term  except  the  first, 
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(C.2.14) 

One  of  equations  C. 2. 9 ,0.2.11  or  C.2.14  is  applied  for  each 
(P,  M,  G)  in  Figure  C.2.1.  The  resulting  equations  are  solved 

simultaneously  to  give  each  of  the  ZC2 .  .  A  much  simpler 

set  of  equations  is  generated  for  Q=S=H=0.  This  condition  applies 
in  the  far  field  case  (see  Appendix  D) .  The  first  of  the  far- 
field  case  equations  is  obtained  by  integrating  the  equation 
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by  parts  with  respect  to  the  dummy  variable  x  to  give 
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Comparing  with  equation  C.2.2,  replacing  N  by  N+2  and  P  by  P-2 
and  separating  all  terms  for  which  i+j+k=2,  the  recursive  relation 
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(C.2.17) 

is  obtained.  This  relation  can  be  applied  for  all  P%1.  2.  second 

far  field  case  equation  is  obtained  by  integrating  equation  C.2.15 
by  parts  with  respect  to  the  dummy  variable  y  to  give 
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(C.2.18) 

Comparing  with  equation  C.2.2,  replacing  N  by  N+2  and  M  by  M-2 
and  separating  all  terms  for  which  i+j+k=2,  the  recursive  rela¬ 
tion 
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This  relation  can  be  applied  for  all  M>1.  A  third  frr  field 
equation 
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is  derived  similarly.  This  relation  can  be  used  for  all  G^l. 

A  constant  P+M+G  plane  is  shown  in  figure  C. 2.2  and  a  strategy 
for  calculating  2c2NP0M0G0  outl^ned* 


-  i 


C- 27 


C.2.20. 

Solve  for  these  six 
points. 

Figure  C.2.2 
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APPENDIX  D 


THE  NUMERICAL  PROCEDURE  FOR  THE  COMPUTATION  OF 


_WW 


_ws 


„sw 


AND  Z 


SS 


Equations  for  the  computation  of  ZQ  ,  z£j  ,  ZQ  and  ZQ 
have  been  developed  in  sections  2.3,  2.4  and  2.5.  These  are 
developed  further  in  this  appendix  to  a  form  that  is  more  con¬ 
venient  for  computational  purposes.  The  partial  sums  are  subject 
to  a  high  degree  of  instability  when  the  distances  are  large 
compared  to  the  sizes  of  the  segments.  For  this  case,  separate 
far  field*  solutions  have  been  developed. 

Equations  2.3.8  and  2.3.12  are  used  for  computing  the 
mutual  impedances  from  Z^.  For  the  purpose  of  computing  Z^, 
it  is  appropriate  to  collect  terms  of  order  (n+p)  together. 
Writing  q  =  n  +p,  and  noting  that 


♦Traditionally,  far  field  implies  a  region  far  enough  so  that  the 
field  on  the  receptor  can  be  considered  a  constant.  This  is  not 
the  implication  here  when  the  term  far  field  is  used.  Instead, 
far  field  implies  a  region  far  enough  so  that  Jl/r ,  the  ratio  of 
segment  size  to  distance  between  the  receiver  and  transmitter 
origins  is  small  enough  that  <£,/r)3  is  negligible  compared  to  1. 
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where  m'  =  m  -  n,  equation  2.3.18  can  be  written  as 
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When  na  «  rno,  all  ZA2v  (s+m. j  (n+a)  ®re  of  the  order  of 
ZA2v(s+m')0  and  the  innermost  summation  results  in  the  sub¬ 
traction  of  almost  equal  quantities.  Therefore,  it  is  appropriate 
and  convenient  to  express  the  quantity 


nf0  (n)  (_1)  ZA2v(s+m*)  (n+a)  ([AU[R])  R0 


(D.  2) 


noting  that  R^  is  of  the  order 


as  a  series  in  R,  and  R„  [noting  that  R.  is  of  the  order 
1  2  l  1 

#  Since  Rq  is  1,  the  last  factor  contributes  nothing 
to  the  expression.  However,  its  introduction  makes  it  easier  to 
generalize  the  equations  being  developed  here  to  wire-to-surf ace 
and  surface-to-surface  equations. 
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The  application  of  the  binomial  theorem  gives 


The  results  of  equation  D.3  are  used  for  rewriting  the  equations 
for  Zq*  in  a  form  that  is  numerically  rugged  when  L<<r.  Thus 
the  quantity  under  consideration  in  equation  D.2  becomes 
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The  quantity  in  the  square  brackets  is  empirically  found  to 
be  zero*  for  80  <  q-  Thus  the  expression  D.4  can  be  written. 


(using  BQ  =  ^j-q)  , 
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A  formal  proof  was  not  attempted  in  this  effort.  However,  that 


the  summation  in  expression  D.2  should  be  of  the  order 
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becomes  obvious  if  one  writes 
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All  factors  except  the  square  brackets  are  of  order  1.  The 
quantity  inside  the  square  parenthesis  is,  by  applying  the 
binomial  theorem  in  reverse. 
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The  quantity  in  the  second  factor  is  of  the  order  1.  Thus  the 
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H  (q+n  o ) 

Z 


nn=0  en=0 


q  +  n  o  ~  e  o 
e  o 


^“-(q  +no  -  eo) 


n  o“2e  o  e  o 

X  R1  R2  Aaq(n -eq)  ZA2v ( s+m* +q+n o ) 0 


(D.  7) 


Thus 


for  <<  r^,  )f  equation  D.l  is  implemented  using 


„WW  2  (-1)°S  exp(- jkr110)  *aS+1  (aJ)^  r“1Q 

Z0  Z i  2j  sin  k  l 

c=l  o 


x  jo  H  *  OiqUl-<-1,ii  jlo 


R1  r110 

q  2 

r 

*.  f 

^(no+q) 

r 

q+no~£o] 

R1  r110 

0 

a 

L 

no=0 

iriioi 

£0=0 

eo  J 

l 

a 

n  o~2e  o 


x  Aaq(no~£o)  ZA2v (s+m' +no+q)  0 


(D.8) 


•  / 


■  •  UM-IP4' 


D-6 


The  above  equation  has  been  written  so  as  to  keep  all  partial 

la 

sums  independent  of  the  ratio  — which  may  vary  considerably 
in  the  computations  depending  upon  the  geometry.  Thus  all 
partial  sums  in  the  above  equation  are  well  behaved. 


Computation  of  wire-to-surf ace  and  surface-to-wire 
impedances  is  performed  using  a  strategy  similar  to  one  used 
for  computation  of  wire-to-wire  impedances.  Equation  2.4.17a 
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It  is  noted  that  the  variable  n  occurs  in  two  of  the  Z__ 
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indices.  This  makes  the  application  of  recursion  relations 
B.2.2,  B.2.3  and  B.2.4  inconvenient  as  the  summation  is  exec¬ 
uted.*  The  summation  is  made  more  convenient  by  writing 
equation  D.9  to  give 
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first  is  to  fill  a  5  dimensional  matrix  with  the  appropriate  values 

of  Zg2  and  compute  according  to  equation  D.9.  This  approach 

is  attractive  but  would  take  approximately  2K  of  core  for  storing 

ZB2  and  approximately  100K  for  storing  Zc2  in  comp- 
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uting  Zq  .  Since  this  is  too  large,  a  second  approach,  that  of 
calculating  the  base  values  of  ZB2  (and  storing  them)  and 

then  using  equations  B.2.2,  B.2.3  and  B.2.4  as  the  summation 
proceeds, must  be  used. 
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or,  writing  q-n+m^n^  as  a  summation  index, 
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As  in  the  wire-to-wire  case,  it  is  appropriate  and  convenient 
to  express  the  quantity 
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in  equation  D.9  as  a  series  in  ( — a — )  when  *xj<<rnin- 
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In  the  wire-to-surface  case,  RQ0  is  1. 

The  last  factor  has  been  included  to  facilitate  generalization 
of  this  development  to  the  surface-to-surface  case.  Expression 
D.ll  can  be  written 
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or  by  a  development  process  similar  to  that  used  in  developing 
(D.  7  )from  (D.  2  )  , 
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It  is  to  be  noted  that  all  partial  suns  in  equation  D.15a  are  veil 
behaved.  Similarly, 
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For  the  surface-to-surface  case,  combining  equations  2.5.6  and 
2.5.9, 
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similar  to  equation  D.10, 
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writing  i1~m1+m2=n2  as  a  summation  index, 
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When  rooo  *s  ^ar(3e  •  tlie  innermost  summations  in  (D.16)  are 
evaluated  by  considering  the  expression 
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This  expression  is  developed  similar  to  the  development  of 
expression  D.7  to  give 
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Thus,  for  the  far -field  case,  equation  D.16  is  written  as 
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MISSION 

of 

Rome  Air  Development  Center 

RAVC  plan*  and  executes  research,  development,  test  and 
selected  acquisition  programs  in  support  o&  Command,  Control 
Communications  and  Intelligence  1C* I)  activities.  Technical 
and  engineeAing  support  utithin  areas  oi  technical  competence 
is  provided  to  ESP  Program  0  ll-ic.es  IPOs)  and  other  ESP 
elements.  The  principal  technical  mission  areas  one 
communications,  electromagnetic  guidance  and  control,  sur¬ 
veillance  oi  ground  and  aerospace  objects,  intelligence  data 
collection  arid  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  sciences,  micromve 
physics  ami  electronic  reliability,  maintainability  and 
compatibility. 
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